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Abstract 

This paper is concerned with the macroscopic behavior of global energy minimizers 
in the three-dimensional sharp interface unscreened Ohta-Kawasaki model of diblock 
copolymer melts. This model is also referred to as the nuclear liquid drop model in 
the studies of the structure of highly compressed nuclear matter found in the crust 
of neutron stars, and, more broadly, is a paradigm for energy-driven pattern forming 
systems in which spatial order arises as a result of the competition of short-range at¬ 
tractive and long-range repulsive forces. Here we investigate the large volume behavior 
of minimizers in the low volume fraction regime, in which one expects the formation of 
a periodic lattice of small droplets of the minority phase in a sea of the majority phase. 
Under periodic boundary conditions, we prove that the considered energy T-converges 
to an energy functional of the limit “homogenized” measure associated with the minor¬ 
ity phase consisting of a local linear term and a non-local quadratic term mediated by 
the Coulomb kernel. As a consequence, asymptotically the mass of the minority phase 
in a minimizer spreads uniformly across the domain. Similarly, the energy spreads 
uniformly across the domain as well, with the limit energy density minimizing the en¬ 
ergy of a single droplets per unit volume. Finally, we prove that in the macroscopic 
limit the connected components of the minimizers have volumes and diameters that 
are bounded above and below by universal constants, and that most of them converge 
to the minimizers of the energy divided by volume for the whole space problem. 
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1 Introduction 


The liquid drop model of the atomic nucleus, introduced by Gamow in 1928, is a classical 
example of a model that gives rise to a geometric variational problem characterized by a 
competition of short-range attractive and long-range repulsive forces Qia (for more recent 
studies, see e.g. for a recent non-technical overview of nuclear models, see, e.g., jlO]). 

In a nucleus, different nucleons attract each other via the short-range nuclear force, which, 
however, is counteracted by the long-range Goulomb repulsion of the constitutive protons. 
Within the liquid drop model, the effect of the short-range attractive forces is captured by 
postulating that the nucleons form an incompressible fluid with fixed nuclear density and 
by penalizing the interface between the nuclear fluid and vacuum via an effective surface 
tension. The effect of Goulomb repulsion is captured by treating the nuclear charge as 
uniformly spread throughout the nucleus. A competition of the cohesive forces which try 
to minimize the interfacial area of the nucleus and the repulsive Goulomb forces that try to 
spread the charges apart makes the nucleus unstable at sufficiently large atomic numbers, 
resulting in nuclear fission |4 11-1^. 
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Figure 1: Nuclear pasta phases in a relativistic mean-field model of low density nuclear 
matter. The panels show a progression from “meatball” (a) to “spaghetti” (b) to “lasagna” 
(c) to “macaroni” (d) to “swiss cheese” (e) phases, which are the numerically obtained 
candidates for the ground state at different nuclear densities. Reproduced from Ref. |25]. 


It is worth noting that the liquid drop model is also applicable to systems of many 
strongly interacting nuclei. Such a situation arises in the case of matter at very high den¬ 
sities, occurring, for example, in the core of a white dwarf star or in the crust of a neutron 
star, where large numbers of nucleons are conhned to relatively small regions of space 
by gravitational forces |14f[I^ . As was pointed out independently by Kirzhnits, Abrikosov 
and Salpeter, at sufficiently low temperatures and not too high densities compressed matter 
should exhibit crystallization of nuclei into a body-centered cubic crystal in a sea of de¬ 
localized degenerate electrons |17f[I^ . At yet higher densities, more exotic nuclear “pasta 
phases” are expected to appear as a consequence of the effect of “neutron drip” 16,20-26 


(for an illustration, see Fig. [^. In all cases, the ground state of nuclear matter is deter¬ 
mined by minimizing the appropriate (free) energy per unit volume of one of the phases 
that contains contributions from the interface area and the Coulomb energy of the nuclei. 

Within the liquid drop model, the simplest way to introduce conhnement is to restrict 
the nuclear fluid to a bounded domain and impose a particular choice of boundary conditions 
for the Coulombic potential. Then, after a suitable non-dimensionalization the energy takes 
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the form 


E{u) := I \Vu\dx + - I I G{x,y){u{x) — u){u{y) — u) dx dy. (1-1) 

Jn 2 Jq 

Here, H C is the spatial domain (bounded), u G 1}) is the characteristic 

function of the region occupied by the nuclear fluid (nuclear fluid density), u G (0,1) is 
the neutralizing uniform background density of electrons, and G is the Green’s function 
of the Laplacian which, in the case of Neumann boundary conditions for the electrostatic 
potential solves 

-A3;G(x,y) = h(x - y) - 1^, (1.2) 

where (5(x) is the Dirac delta-function. The nuclear fluid density must also satisfy the global 
electroneutrality constraint: 


7^ j udx = u. 

I“l Jo. 


(1.3) 


In writing (1.1) we took into account that because of the scaling properties of the Green’s 


function one can eliminate all the physical constants appearing in (1.1) by choosing the 
appropriate energy and length scales. 


It is notable that the model in (1.1)-(1.3) also appears in a completely different physical 


context, namely, in the studies of mesoscopic phases of diblock copolymer melts, where it 
is referred to as the Ohta-Kawasaki model 27-29 . This is, of course, not surprising. 


considering the fundamental nature of Goulomb forces. In fact, the range of applications 


of the energy in (1.1) goes far beyond the systems mentioned above (for an overview, 
see and references therein). Importantly, the model in (1.1) is a paradigm for the 
energy-driven pattern forming systems in which spatial patterns (global or local energy 
minimizers) form as a result of the competition of short-range attractive and long-range 
repulsive forces. This is why this model and its generalizations attracted considerable 
attention of mathematicians in recent years (see, e.g., |31f|49| , this list is certainly not 
exhaustive). In particular, the volume-constrained global minimization problem for (1.1) 


in the whole space with no neutralizing background, which we will also refer to as the 
“self-energy problem”, has been investigated in p4 37 45 49 . 


A question of particular physical interest is how the ground states of the energy in (1.1) 


behave as the domain size tends to inhnity. In [^, Alberti, Ghoksi and Otto showed that 
in this so-called “macroscopic” limit the energy becomes distributed uniformly throughout 
the domain. Another asymptotic regime, corresponding to the onset of non-trivial mini¬ 
mizers in the two-dimensional screened version of (1.1) was studied in (33], where it was 


shown that at appropriately low densities every non-trivial minimizer is given by the char¬ 
acteristic function of a collection of nearly perfect, identical, well separated small disks 
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(droplets) uniformly distributed throughout the domain (see also |38] for a related study 
of almost minimizers). Further results about the fine properties of the minimizers were 
obtained via two-scale F-expansion in (^, using the approach developed for the studies of 
magnetic Ginzburg-Landau vortices |50| (more recently, the latter was also applied to three- 
dimensional Coulomb gases |51j). In particular, the method of |39] allows, in principle, to 
determine the asymptotic spatial arrangement of the droplets of the low density phase via 
the solution of a minimization problem involving point charges in the plane. It is widely 
believed that the solution of this problem should be given by a hexagonal lattice, which in 
the context of type-II superconductors is called the “Abrikosov lattice” |52| . Proving this 
result rigorously is a formidable task, and to date such a result has been obtained only 
within a much reduced class of Bravais lattices 1^,53 . 


It is natural to ask what happens with the low density ground state of the energy in (1.1) 
as the size of the domain goes to infinity in three space dimensions. As can be seen from 
the above discussion, the answer to this question bears immediate relevance to the structure 
of nuclear matter under the conditions realized in the outer crust of neutron stars. This is 
the question that we address in the present paper. On physical grounds, it is expected that 
at low densities the ground state of such systems is given by the characteristic function of 
a union of nearly perfect small balls (nuclei) arranged into a body-centered cubic lattice 
(known to minimize the Coulomb energy of point charges among body-centered cubic, face- 
centered cubic and hexagonal close-packed lattices |54f|56| ). The volume of each nucleus 
should maximize the binding energy per nucleon, which then yields the nucleus of an isotope 
of nickel. 

Our results concerning the minimizers of (1.1) proceed in that direction, but are still 
far from rigorously establishing such a detailed physical picture. One major difficulty has 
to do with the lack of the complete solution of the self-energy problem [^[48] . Assuming 
the solution of this problem, whenever it exists, is a spherical droplet, a mathematical 
conjecture formulated by Choksi and Peletier |35| and a universally accepted hypothesis 
in nuclear physics, we indeed recover spherical nuclei whose volume minimizes the self¬ 
energy per unit nuclear volume (which is equivalent to maximizing the binding energy per 
nucleon in the nuclear context). The question of spatial arrangement of the nuclei is another 
major difficulty related to establishing periodicity of ground states of systems of interacting 
particles, which goes far beyond the scope of the present paper. Nevertheless, knowing that 


the optimal droplets are spherical should make it possible to apply the techniques of |50 51 


to relate the spatial arrangement of droplets to that of the minimizers of the renormalized 
Coulomb energy. 

In the absence of the complete solution of the self-energy problem, we can still establish, 
although in a somewhat implicit manner, the limit behavior of the minimizers of (1.1)-(1.31 
in the case Q = T^, where is the three-dimensional torus with sidelength i, as i ^ oo, 
provided that u also goes simultaneously to zero with an appropriate rate (low-density 
regime). We do so by establishing the F-limit of the energy in (1.1), with the notion of 
convergence given by weak convergence of measures (for a closely related study, see |38] ). 


5 











The limit energy is given by the sum of a constant term proportional to the volume occupied 
by the minority phase (which is also referred to as “mass” throughout the paper) and the 
Coulombic energy of the limit measure, with the proportionality constant in the first term 
given by the minimal self-energy per unit mass among all masses for which the minimum of 
the self-energy is attained. Importantly, the minimizer of the limit energy (which is strictly 
convex) is given uniquely by the uniform measure. Thus, we establish that for a minimizer 


of (1.1)-(1.3) the mass in the minority phase spreads (in a coarse-grained sense) uniformly 
throughout the spatial domain and that the minimal energy is proportional to the mass, 
with the proportionality constant given by the minimal self-energy per unit mass (compare 
to ( 32 ]). We also establish that almost all the “droplets”, i.e., the connected components of 
the support of a particular minimizer, are close to the minimizers of the self-energy with 
mass that minimizes the self-energy per unit mass. 

Mathematically, it would be natural to try to extend our results in two directions. The 


first direction is to consider exactly the same energy as in (1.1) in higher space dimensions. 


Here, however, we encounter a difficulty that it is not known that the minimizers of the 
self-energy do not exist for large enough masses. Such a result is only available in three 
space dimensions for the Coulombic kernel 37,45 . In the absence of such a non-existence 


result one may not exclude a possibility of a network-like structure in the macroscopic limit. 


Another direction is to replace the Coulombic kernel in (1.1) with the one corresponding 
to a more general negative Sobolev norm. Here we would expect our results to still hold 
in two space dimensions. Furthermore, the physical picture of identical radial droplets 
in the limit is expected for sufficiently long-ranged kernels, i.e., those kernels that satisfy 
G{x,y) ~ |x — for |x — 7/| <C 1, with 0 < a <C 1 36,^. Note that although a similar 


characterization of the minimizers for long-ranged kernels exists in higher dimensions as 
well |^|46| , these results are still not sufficient to be used to characterize the limit droplets, 
since they do not give an explicit interval of existence of the minimizers of the self-interaction 
problem. Also, since the non-existence result for the self-energy with such kernels is available 
only for a < 2 37 , our results may not extend to the case of a > 2 in dimensions three 


and above. 

Finally, a question of both physical and mathematical interest is what happens with 
the above picture when the Coulomb potential is screened (e.g., by the background density 
fluctuations). In the simplest case, one would replace (1.2) with the following equation 
defining G: 


-A^G{x, y) + k^G{x, y) = 6(x - y), 


(1.4) 


where k > 0 is the inverse screening length, and the charge neutrality constraint from (1.3) 


is relaxed. Here a bifurcation from trivial to non-trivial ground states is expected under 


suitable conditions (in two dimensions, see [M 38 ^). We speculate that in certain limits 
this case may give rise to non-spherical droplets that minimize the self-energy. Indeed, in 
the presence of an exponential cutoff at large distances, it may no longer be advantageous 
to split large droplets into smaller disconnected pieces, and the self-energy minimizers for 
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arbitrarily large masses may exist and resemble a “kebab on a skewer”. In contrast to the 
bare Coulomb case, in the screened case the energy of such a kebab-shaped conhguration 
scales linearly with mass. Note that this conhguration is reminiscent of the pearl-necklace 


morphology exhibited by long polyelectrolyte molecules in poor solvents 157,58 


Organization of the paper. In Sec. we introduce the specihc model, the scaling 
regime considered, the functional setting and the heuristics. In this section, we also discuss 
the self-energy problem and mention a result about attainment of the optimal self-energy 
per unit mass. In Sec. we hrst state a basic existence and regularity result for the 
minimizers (Theorem |3.1[) and give a characterization of the minimizers of the whole space 


problem that also minimize the self-energy per unit mass (Theorem 3.2). We then state 


our main T-convergence result in Theorem 3.3 In the same section, we also state the 


consequences of Theorem 3.3 to the asymptotic behavior of the minimizers in Corollary 


3.4, as well as Theorem 3.5 about the uniform distribution of energy in the minimizers 


and Theorem |3.6| that establishes the multidroplet character of the minimizers. Section]^ 
is devoted to generalized minimizers of the self-energy problem, where, in particular, we 


obtain existence and uniform regularity for minimizers in Theorem 4.5 and Theorem |4.7| 
This section also establishes a connection to the minimizers of the whole space problem with 
a truncated Coulombic kernel and ends with a characterization of the optimal self-energy 
per unit mass in Theorem 4.15 Section contains the proof of the T-convergence result 
of Theorem 3.3 and of the equidistribution result of Theorem |3.5| Section establishes 
uniform estimates for the problem on the rescaled torus, where, in particular, uniform 
estimates for the potential are obtained in Theorem 6.9 Section presents the proof of 
Theorem |3.6[ Finally, some technical results concerning the limit measures appearing in 
the T-limit are collected in the Appendix. 


Notation. Throughout the paper , BV, A4 denote the usual spaces 

of Sobolev functions, functions of bounded variation, Lebesgue functions, functions with 
continuous derivatives up to order k, compactly supported functions with continuous deriva¬ 
tives up to order k, functions with Holder-continuous derivatives up to order k for a G (0,1), 
and the space of hnite signed Radon measures, respectively. We will use the symbol |Vu| 
to denote the Radon measure associated with the distributional gradient of a function of 
bounded variation. With a slight abuse of notation, we will also identify Radon measures 
with the associated, possibly singular, densities (with respect to the Lebesgue measure) on 
the underlying spatial domain. For example, we will write v = |Vn| and dv{x) = |Vn(x)| dx 
to imply u G A4(H) and i/(H') = |Vu|(H') = |Vu| dx, given u G BV(fi) and C H. 
For a set / C N, denotes cardinality of I. The symbol XF always stands for the char¬ 
acteristic function of the set F, and \F\ denotes its Lebesgue measure. We also use the 
notation (ug) G to denote sequences of functions G Ae as e = £n ^ 0, where Ae are 
admissible classes. 
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2 Mathematical setting and scaling 


Variational problem on the unit torus. Throughout the rest of this paper the spatial 
domain 0 in (1.1) is assumed to be a torus, which allows us to avoid dealing with boundary 
effects and concentrate on the bulk properties of the energy minimizers. We dehne T : = 
to be the flat three-dimensional torus with unit sidelength. For e > 0, which should 


be treated as a small parameter, we introduce the following energy functional: 


E£{u) := e f \Vu\dx + - f {u — Ue){—A) ^{u — Ue)dx, (2-1) 

JT 2 Jf 

where the hrst term is understood distributionally and the second term is understood as the 
double integral involving the periodic Green’s function of the Laplacian, with u belonging 
to the admissible class 


where 


A := \u£ BV{T-,{0,1}) 


Udx = Ue 


IT 


( 2 . 2 ) 


u, 


:=A£2/3, 


(2.3) 


with some hxed A > 0. The choice of the scaling of Us with e in (2.3) will be explained 
shortly. To simplify the notation, we suppress the explicit dependence of the admissible 
class on A, which is hxed throughout the paper. 


It is natural to dehne for u € Ae the measure by 

d^i.{x) := e~‘^^^u{x) dx. 


(2.4) 


In particular, is a positive Radon measure and satishes Therefore, on a 

suitable sequence as e —)• 0 the measure /Xg converges weakly in the sense of measures to a 
limit measure /r, which is again a positive Radon measure and satishes ^(T) = A. 


Function spaces for the measure and potential. 

in (2.11 is given by 


In terms of fj,£ the Coulombic term 


1 

2 



n£)(-A) ^{u 


Ue) dx 



y)dne{x)dne{y), 


(2.5) 


where G is the periodic Green’s function of the Laplacian on T, i.e., the unique distributional 
solution of 


—AG(x) = h(x) — 1, 
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G{x) dx = 0. 


( 2 . 6 ) 






If the kernel G in (2.5) were smooth, then one wonld be able to pass directly to the limit 


in the Conlombic term and obtain the corresponding convolntion of the kernel with the 
limit measnre. This is not possible dne to the singnlarity of the kernel at {x = y}. In 
fact, the donble integral involving the limit measnre may be strictly less than the liminf 
of the seqnence, and the defect of the limit is related to a non-trivial contribntion of the 
self-interaction of the connected components of the set {u = 1} and its perimeter to the 
limit energy. 

On the other hand, the singnlar character of the kernel provides control on the regnlarity 
of the limit measnre fi. To see this, we dehne the electrostatic potential G by 


which solves 


:= / G{x -y)dfie{y), 

Jt 


/ Vip -Vvedx = / (pdpe — ^ ipdx 

h Jt Jt 


\fp E C~(T). 


By (2.4), we can rewrite the corresponding term in the Conlombic energy as 


G{x - y) dpe{x) dpeiv) = / Vedpe= / \^Ve\ dx. 

IT JT Jt Jt 


(2.7) 


( 2 . 8 ) 


(2.9) 


Hence, if the left-hand side of (2.9) remains bonnded as e —)> 0, and since dx = 0, the 
seqnence is nniformly bonnded in 77^ (T) and hence weakly convergent in H^{T) on a 
snbseqnence. 


By the above discnssion, the natnral space for the potential is the space 


H := E 77^(T) : J ndx = o| with 


Mn ■= 


|Vnp dx 


1/2 


( 2 . 10 ) 


The space T-L is a Hilbert space together with the inner prodnct 


{u,v)'n:= / Vu-Vvdx 
Jt 


Mu, V eTi. 


( 2 . 11 ) 


The natnral class for measnres pe to consider is the class of positive Radon measnres on T 
which are also in Ti', the dnal of %. More precisely, let A1“'"(T) C A4(T) be the set of all 
positive Radon measnres on T. We dehne the snbset A4'''(T) n %' of A4+(T) by 

M + (T)n'H' = {/xEM+(T) : jpdyi<G\\p\\H E H n G(T)|, (2.12) 
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for some C > 0. This is the set of positive Radon measures which can be understood as 
continuous linear functionals on %. Note that /r G A^'’'(T) satisfies // G A^'*'(T) r\l-L' if and 
only if it has hnite Coulombic energy, i.e. 


T JT 


G{x — y) dfi{x) dn{y) < oo, 


(2.13) 


with the convention that G'(O) = + 00 . The proof of this characterization and related facts 
about n %' are given in the Appendix. 


The whole space problem. We will also consider the following related problem, 
formulated on M^. We consider the energy 

Eoo{u) := f \Vu\dx + -^ f f dxdy. (2-14) 

Jr3 Svr J^3 J^3 \x — y\ 

The appropriate admissible class for the energy E^o in the present context is that of con- 
hgurations with prescribed “mass” m > 0: 

:= {0,1}) : J udx = m^. (2.15) 

For a given mass m > 0, we dehne the minimal energy by 

e{m) := inf Eao{u). (2.16) 

ueAao(m) 


The set of masses for which the inhmum of E^o in Aoo{m) is attained is denoted by 

X := |m > 0 : 3 G ^oo(?7r), .Eoo(um) = e(m)| , (2-17) 

The minimization problem associated with (|2.14 ) and (2.15) was recently studied by two 
of the authors in |37]. In particular, by 37 Theorem 3.3] the set X is bounded, and 


by |37[ Theorems 3.1 and 3.2] the set X is non-empty and contains an interval around the 
origin. 

For m > 0, we also dehne the quantity (with the convention that /(O) ;= -|-oo) 

e{m) 


f{m) : = 


m 


(2.18) 


which represents the minimal energy for (2.14) and (2.15) per unit mass. By 37 Theorem 
3.2] there is a universal fho > 0 such that f{m) is obtained by evaluating E^o on a ball of 
mass m for all m < friQ. After a simple computation, this yields 

/(m) = -h 3^/^ • 2"^/^ • 10"^ • for all 0 < m < mo- (2.19) 
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Note that obviously this expression also gives an a priori upper bound for f{m) for all 
m > 0. In addition, by [37{ Theorem 3.4] there exist universal constants C,c> 0 such that 

( 2 . 20 ) 


c < f{m) < C for all m > mo¬ 
lt was conjectured in that X = [0,mo] and that mg = urd, where 

dOvr 


rrici ■■= 


(^2^/3 + 2-1/3 _ ~ 44.134. 


( 2 . 21 ) 


The quantity md is the maximum value of m for which a ball of mass m has less energy 
than twice the energy of a ball with mass ^m. However, such a result is not available at 
present and remains an important challenge for the considered class of variational problems 


(for several related results see [M 44 47]). 
Finally, we dehne 

f* := inf /(m) and 

m£X 


X* :={m* GX : /(m*) =/*} 


( 2 . 22 ) 


Observe that in view of (2.19) and (2.20) we have f* G (0,oo). Also, as we will show in 
Theorem 3.2 the set X* is non-empty, i.e., the minimum of /(m) over X is attained. In 


fact, the minimum of /(m) over X is also the minimum over all m G (0, oo) (see Theorem 


4.15). Note tha t this result was also independently obtained by Frank and Lieb in their 
The set X* of masses that minimize the energy E^o per unit mass and 


49 


recent work 

the associated minimizers (which in general may not be unique) will play a key role in the 


analysis of the limit behavior of the minimizers of E^. Note that if f{m) were given by (2.19) 
and X = [0, md], then we would have explicitly X* = {lOvr} and f* = 3^/3 • 2-^/3 • 5 - 1/3 ~ 


2.29893. On the other hand, in view of the statement following (2.19), this value provides 
an a priori upper bound on the optimal energy density. 


Macroscopic limit & heuristics. The limit e —)• 0 with A > 0 hxed is equivalent to the 
limit of the energy in (1.1) with Q = T^, where := is the torus with sidelength 

£ > 0, as £ —^ 00 . Indeed, introducing the notation 


1 




/( 

JTe 


Ei{u) := / \Vu\dx + - / {u-ue){-A) \u-ui)dx, 


for the energy in (1.1) with H = T^, and taking U£ = X£ ^ and h G where 


Ai := u G i?H(T£; {0,1}) : / udx = X£ 

Jjt 


(2.23) 


(2.24) 


it is easy to see that u{x) := u{£x) belongs to Ae with Us = Ae^/3 for e = £ ^, and we have 
E^{u) = £~^Ei(u). It then follows that the two limits 00 and e —>■ 0 are equivalent. 


Note that the full space energy E^c is the formal limit of (2.23) for £ —)> 00 . 
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The choice of the scaling of u,. with e is determined by the balance of far-held and near- 
held contributions of the Coulomb energy. Heuristically, one would expect the minimizers 


of the energy in (2.1) to be given by the characteristic function of a set that consists of 
“droplets” of size of order i? <C 1 separated by distance of order d satisfying R d 1 
(for evidence based on recent molecular dynamics simulations, see also |26]). Assuming 


that the volume of each droplet scales as (think, for example, of all the droplets as non¬ 


overlapping balls of equal radius and with centers forming a periodic lattice), from (2.23) 
we hnd for the surface energy, self-energy and interaction energy, respectively, for a single 
droplet: 

A'surf ^ £R , T'self ~ R i 

Equating these three quantities and recalling that R?/d? ~ hg, we obtain 

R ~ d ~ hj 


R^ 


Eint ~ 

(2.25) 

^ ^ Ue, we obtain 


r.2/3 

(2.26) 


which leads to (2.3). Note that, in some sense, this is the most interesting low volume 
fraction regime that leads to inhnitely many droplets in the limit as e —>■ 0, since both the 
self-energy of each droplet and the interaction energy between different droplets contribute 
comparably to the energy. For other scalings one would expect only one of these two terms 
to contribute in the limit, which would, however, result in loss of control on either the 
perimeter term or the Coulomb term as e —)> 0 and, as a consequence, a possible change in 
behavior. Let us note that a different scaling regime, in which Ue = 0(e), leads instead to 
finitely many droplets that concentrate on points as e —)■ 0 p4] , while for Ug = 0(1) one 
expects phases of reduced dimensionality, such as rods and slabs (see Fig. [^. 


3 Statement of the main results 


We now turn to stating the main results of this paper concerning the asymptotic behavior of 


the minimizers or the low energy configuration of the energy in (2.1) within the admissible 


class in (2.2). Existence of these minimizers is guaranteed by the following theorem. 


Theorem 3.1 (Minimizers: existence and regularity). For every A > 0 and every 0 < 
e < there exists a minimizer rtg G Ag of given by (2.1) with Us given by (2.3). 

Furthermore, after a possible modification of Us on a set of zero Lebesgue measure the 
support of Ue has boundary of class C°°. 


Proof. The proof of Theorem |3.1| is fairly standard. We present a few details below for the 
sake of completeness. 

By the direct method of the calculus of variations, minimizers of the considered problem 
exist for all e > 0 as soon as the admissible class As is non-empty, in view of the fact 
that the first term is coercive and lower semicontinuous in BV{T), and that the second 
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term is continuous with respect to the L^(T) convergence of characteristic functions. The 
admissible class is non-empty if and only if e < 

Holder regularity of minimizers was proved in |59[ Proposition 2.1], where it was shown 
that the essential support of minimizers has boundary of class Smoothness of the 


boundary was established in |43[ Proposition 2.2] (see also the proof of Lemma 4.4 below 
for a brief outline of the argument in a closely related context). □ 

In view of the regularity statement above, throughout the rest of the paper we always 
choose the regular representative of a minimizer. 


We proceed by giving a characterization of the quantity f* defined in (2.22) as the 
minimal self-energy of a single droplet per unit mass, i.e., as the minimum of f{m) over I. 

Theorem 3.2 (Self-energy: attainment of optimal energy per unit mass). Let f* be defined 
as in (2.22). Then there exists m* G X such that f* = f{m*). 

With the result in Theorem 


3.2[ we are now in the position to state our main result on the 
T-limit of the energy in (2.1), which can be viewed as a generalization of 38, Theorem 1]. 


Theorem 3.3 (P-convergence). For a given A > 0, let he defined by (2.1) with Ue given 

0 we have ^ 


by (2.3). Then as e 


Eoifi) := Xr + ^ 


Eq, where 

G{x -y)dn{x)dix{y), 


(3.1) 


T JT 


and fi G A4+(T) n TL' satisfies /u(T) = A. More precisely, 

i) (Compactness and T-liminf inequality) Let (up) G As be such that 

limsupe“^/^ille(ne) < oo, 

£—^0 


(3.2) 


and let p,s and Vs be defined in (2.4) and (2.7), respectively. Then, upon extraction of 
a subsequence, we have 


/Te ^ /r m A4(T), Vs ^ v in TL, (3.3) 

as e —>■ 0, for some [i G {T)r\TL' with ^(T) = A, the function v has a representative 
in L^{T,dp,) given by 

v{x)= [ G{x -y)dn{y), (3.4) 

Jt 

and 

ImiuiiEsius) > E^iy). (3.5) 
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ii) (V-limsup inequality) For any measure n G A^'*“(T) r\%' with /u(T) = A there exists a 
sequence (us) G As such that (3.3) and (3.4) hold as e — )• 0 for pLs CL^d Ve defined in 


(2.4) and (2.7), and 


lim sup e 

£—^0 




(3.6) 


Note that the weak convergence of measures was recently identihed in p8) (see also 
as a suitable notion of convergence for the studies of the F-limit of the two-dimensional 


version of the energy in (2.1) 


Observe also that the limit energy Eq is a strictly convex functional of the limit measure 
and, hence, attains a unique global minimum. By direct inspection, Eq is minimized by 
H = fiQ, where d/xo := Xdx. Thus, the quantity f* plays the role of the optimal energy 
density in the limit e —)■ 0. 

The remaining results are concerned with sequences of minimizers. We will hence assume 
that the functions (us) G As are minimizers of the functional Es- In this case, we can give 
a more precise characterization for the asymptotic behavior of the sequence. We hrst note 


the following immediate consequence of Theorem 3.3 for the convergence of sequences of 
minimizers. 


Corollary 3.4 (Minimizers: uniform distribution of mass). For A > 0, let (us) G As be 


minimizers of E;,, and let and Vs be defined in (2.4) and (2.7), respectively. Then 


Te ho in A4(T), Vs ^ 0 in Fi, 


(3.7) 


where dpLQ = Xdx, and 


e-^/^Efius) ^ Xf, 


(3.8) 


where f* is as in (2.22), as e —)> 0. 


The formula in (3.8) suggests that in the limit the energy of the minimizers is dominated 


by the self-energy, which is captured by the minimization problem associated with the 
energy E^o dehned in (2.14). Therefore, it would be natural to expect that asymptotically 
every connected component of a minimizer is close to a minimizer of ifoo under the mass 
constraint associated with that connected component. Note that in a closely related problem 


in two space dimensions such a result was established in l33l for minimizers, and in 138,39 


for almost minimizers. The situation is, however, unique in two space dimensions, because 
the non-local term in some sense decouples from the perimeter term. Hence, the minimizers 
behave as almost minimizers of the perimeter and, therefore, are close to balls. In three 
dimensions, however, the perimeter and the non-local term of the self-energy Eoo are fully 
coupled, and, therefore, rigidity estimates for the perimeter functional alone may not 
be sufficient to conclude about the “shape” of the minimizers. Nevertheless, we are able to 
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prove a result about the uniform distribution of the energy density of the minimizers as 
e —7- 0 in the spirit of that of |32] . For a minimizer Ue, the energy density is associated with 
the Radon measure dehned by 


dve ■= £ + ^e^^^UeV^dx, 


(3.9) 


where is given by (2.7) and (2.4). Furthermore, we are able to identify the leading order 
constant in the asymptotic behavior of the energy density. 


Theorem 3.5 (Minimizers: uniform distribution of energy). For A > 0, let (ug) G Ae 
minimizers of and let he defined in (3.9). Then 


be 


Ve uq in A4(T) as e —)■ 0, 


(3.10) 


where duQ = Xf*dx and f* is as in (2.22). 


Finally, we characterize the connected components of the support of the minimizers of 
Ei; and show that almost all of them approach, on a suitable sequence as e —>• 0 and after 
a suitable rescaling and translation, a minimizer of Floo with mass in the set I*. 


Theorem 3.6 (Minimizers: droplet structure). For A > 0, let (ue) G Ag be regular rep¬ 
resentatives of minimizers of Eg, let Ng be the number of the connected components of the 
support of Ug, let Ug^k £ {0,1}) be the characteristic function of the k-th connected 

component of the support of the periodic extension of Ug to the whole of modulo trans¬ 
lations in and let Xg^k £ supp(n£^fc). Then there exists Eq > 0 such that the following 
properties hold: 


i) There exist universal constants C,c> 0 such that for all e < Eq we have 

II^£||l°°(t) < C / Ug^kdx>c£, 


(3.11) 


where Vg is given by (2.7). 


a) There exist universal constants C,c > 0 such that for all e < eq we have 
supp(ue^fc) ^ Bfj^i/sixg^k) and cXe~^^^ < Ng < CXe~^^^. 


(3.12) 


Hi) There exists Ng < Ng with Ng/Ng —)• 1 as e —>■ 0 and a subsequence —>• 0 such 
that for every kn < Ng^ the following holds: After possibly relabeling the connected 
components, we have 

Un^u m L^(M^), (3.13) 


1 /3 

where Un{x) := Ug^^k„{£n (x + Xg^^k„)), and u is a minimizer of E^a over Aocim,*) 
for some m* G T*, where T* is defined in (2.22). 
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The significance of this theorem lies in the fact that it shows that all the connected 
components of the support of a minimizer for sufficiently small e look like a collection of 
droplets of size of order separated by distances of order on average. In particular, 
the conclusion of the theorem excludes configurations that span the entire length of the 
torus, such as the “spaghetti” or “lasagna” phases of nuclear pasta (see Fig. Q. Thus, 
the ground state for small enough e > 0 is a multi-droplet pattern (a “meatball” phase). 
Furthermore, after a rescaling most of these droplets converge to minimizers of the non-local 
isoperimetric problem associated with E^o that minimize the self-energy per unit mass. 


4 The problem in the whole space 


In this section, we derive some results about the single droplet problem from (2.14)-(2.15) 
and the rescaled problem from (2.23)-(2.24). 


4.1 The truncated energy 

For reasons that will become apparent shortly, it is helpful to consider the energies where 
the range of the nonlocal interaction is truncated at certain length scale R. We choose a 
cut-off function r] G with ri'{t) < 0 for all t G M, r]{t) = 1 for all t < 1 and r]{t) = 0 

for all t >2. In the following, the choice of r] is fixed once and for all, and the dependence 
of constants on this choice is suppressed to avoid clutter in the presentation. For i? > 0, 
we then define rjn G by r]ji(x) := r]{\x\/R). For u G we consider the 

truncated energy 


:= [ \Vu\dx + 
4r3 


Vr{x - y)u{x)u{y) 
8vr|x - y\ 


dx dy. 


(4.1) 


This functional will be useful in the analysis of the variational problems associated with 
Eoo and E^. We recall that by the results of 


61 


for each ii > 0 and each m > 0 there 
exists a minimizer of E^ in Aao{iR)- Furthermore, after a possible redefinition on a set 
of Lebesgue measure zero, its support has boundary of class and consists of finitely 

many connected components. Below we always deal with the representatives of minimizers 
that are regular. 


The following uniform density bound for minimizers of the energy is an adaption of (37{ 
Lemma 4.3] for the truncated energy E^ and generalizes the corresponding bound for 
minimizers of E^o- 


Lemma 4.1 (Density bound). There exists a universal constant c > 0 such that for every 
minimizer u G Aoo{m) of E^ for some R,m > 0 and any xq G T we have 


I Br{xo) 


udx> cr^ for all r < min(l, m^'^^). 


(4.2) 
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Proof. The claim follows by an adaption of the proof of (37| Lemma 4.3] to our truncated 


energy E^. Indeed, it is enough to show that the statement of 37, Lemma 4.2] holds with 


Eao replaced by E^. The proof of this statement needs to be modihed, since the kernel 
in the dehnition of E^ is not scale-invariant. We sketch the necessary changes, using the 
same notation as in 1371. 


The construction of the sets E and F proceeds as in the proof of 37 Lemma 4.3]. 


The upper bound 37 Eq. (4.6)] still holds since E^{u) < Eao{u). Related to the cut-off 


function in the definition of E^, we get an additional term in the right-hand side of the 
hrst line of |37[ Eq. (4.6)], which is of the form 


EFi JiFj_ 


vr{x - y) 

\x - y|" 


dx dy — i 


2n—a 


r]Rix - y) 


dx dy 


JfJfi k-yh 

^ 2 n-a f f -y)- mix - y) 

Jfi Jfi 


|x - yf 


dx dy < 0, 


(4.3) 


since i > 1 and since the function y is monotonically decreasing (note that a = 1 in 
our case). Since this term has a negative sign, [37] Eq. (4.6)] still holds. The rest of the 
argument then carries through unchanged. □ 


The following lemma establishes a uniform diameter bound for the minimizers of E^. 
The idea of the proof is similar to the one in 45 Lemma 5]. 


Lemma 4.2 (Diameter bound). There exist universal eonstants Rq > 0 and Dq > 0 sueh 
that for any R> Rq, any m > 0 and for any minimizer u G AooifR) of E^, the diameter 
of each connected component Eq o/supp(rt) is bounded above by Dq. 


Proof. Let Eq be a connected component of the support of u with tuq := |To|. Since u is 
a minimizer, xFq is also a minimizer of E^over Aoc{fno)- Indeed, if not, replacing u with 
R~XFo where Xp^ is a minimizer of E^ over AooimQ) translated sufficiently far from 

the support of u would lower the energy, contradicting the minimizing property of u. 


We may assume without loss of generality that R > 2 and diam Eq > 2. Then there 
is G N such that 2N < diam Eq < 2{N + 1). In particular there exist xq, ... ,xj\f G Eq 
such that \xk — xq\ = 2k for every 1 < k < N and, therefore, the balls Bi(xk) are mutually 
disjoint. If mo < 1, then by Lemma 


some universal c > 0. Therefore, 


4.1 


1 /3 

we have |To H Br{xk)\ > cmo for r = mg' < 1 and 


N 

mo>'^ |Fo n Br{xk)\ > citiqN, 
k=l 


(4.4) 


implying that N < Nq for some universal Nq > 1 and, hence, diam Tq < 2(A^o + !)■ If) on 
the other hand, tuq > 1, then by Lemma 4.1 we have \Eq n Bi(xk)\ > c for some universal 
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c > 0. By monotonicity of the kernel in |x — y|, we get 

r I + >Cmin{logAr.logJi;}. 

7F,nB.(«)7F„\B.(.„) If- 91 2t + 2 - 

for some universal C > 0. Hence, if R and N are sufficiently large, then it is energetically 
preferable to move the charge in Bi{xq) sufficiently far from the remaining charge. More 
precisely, consider u = u — XFonBi{xo)~^XFor\Bi{xo){'~^b), for some 6 G with \b\ sufficiently 
large. Then u G .4,oo(?7ro) and 

E^{u) < E^{u) + 4:F - iC'min{logA^,log/?} < 0, (4.5) 

for all R > Rq and N > Nq for some universal constants Rq >2 and Nq > 1. Therefore, 
minimality of u implies that N < Nq whenever R> Rq and hence diam Fq < 2(A^o + l)- ^ 


4.2 Generalized minimizers of 

We begin our analysis of Sqo by introducing the notion of generalized minimizers of the 
non-local isoperimetric problem. 

Definition 4.3 (Generalized minimizers). Given m > 0, we eall a generalized minimizer 
of Eao in Aoo{m) a eollection of functions (ui,... ,ui\f) for some N such that Ui is a 
minimizer of E^o over AooiBXi) with uii = fjj, Ui dx for all i G {1,... ,N}, and 

N N 

171 = '^^ mi and e(m) = e(mi). (4.6) 

i=l i=l 


Clearly, every minimizer of E^ in Aoo{m) is also a generalized minimizer (with = 1). As 
was shown in 


37 


however, minimizers of E^o in Aocina) may not exist for a given m > 0 
because of the possibility of splitting their support into several connected components and 
moving those components far apart. As we will show below, this possible loss of compactness 
of minimizing sequences can be compensated by considering characteristic functions of sets 
whose connected components are “infinitely far apart” and among which the minimum of 
the energy is attained (by a generalized minimizer with some > 1). We also remark that, 
if (ui,..., uat) is a generalized minimizer, then, as can be readily seen from the definition, 
any sub-collection of ufs is also a generalized minimizer with the mass equal to the sum of 
the masses of its components. 


We now proceed to demonstrating existence of generalized minimizers of E^o for all 
m > 0. We start by stating the basic regularity properties of the minimizers of E^o and 
the associated Euler-Lagrange equation. 
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Lemma 4.4 (Regularity and Euler-Lagrange equation). For m > 0, let u be a minimizer 
of Eoo in Aoo{ni), and let F = supp(m). Then, up to a set of Lebesgue measure zero, the 
set F is a bounded connected set with boundary of class C°°, and we have 


2k{x) + vf{x) = Xp for x G dF, 


(4.7) 


where Xp £ M is a Lagrange multiplier, k{x) is the mean curvature of dF at x (positive if 
F is convex), and 


vp{x) 


gL [ 

dvr Jp \x-y\' 


(4.8) 


Moreover, if m £ [mQ,mi] for some 0 < mo < mi, then vp £ (^^’“(M^) and dF is of class 
for all a £ (0,1), uniformly in m. 


Proof. From Proposition 2.1 and Lemma 4.1] it follows that, up to a set of Lebesgue 
measure zero, the set F is bounded and connected, and dF is of class ^ Since the 

function vp is the unique solution of the elliptic problem —An = XF with v{x) -£ 0 for 
|x| —7- oo, by Lemma 4.4] and elliptic regularity theory 62 it follows that vp £ (^^’"(M^) 
for all a £ (0,1), uniformly in m G [mo,mi]. The Euler-Lagrange equation (4.7) can 
be obtained as in |63[ Theorem 2.3] (see also (^[^ ). Further regularity of dF follows 
from Proposition 2.1] and |43[ Proposition 2.2]. □ 

Similarly, if (tti,..., uat) is a generalized minimizer of E^o and T) := supp (ttj) for 
iG{l,...,A^}, the following Euler-Lagrange equation holds: 


2Ki(x) + 7 *- / I I = A x£ dFi, 
4vr Jf, \x - y\ 


(4.9) 


where Ki is the mean curvature of dEi (positive if T) is convex) and A G M is a Lagrange 
multiplier independent of i. 

In contrast to minimizers, generalized minimizers of ifoo in Aloo(ni) exist for all m > 0: 


Theorem 4.5 (Existence of generalized minimizers). For any m £ (0,oo) there exists a 
generalized minimizer (ui,... ,U]\f) of E^o in AoD{m). Moreover, after a possible modifi¬ 
cation on a set of Lebesgue measure zero, the support of each component Ui is bounded, 
connected and has boundary of class C°°. 


Proof. We may assume that m > mo, where mo > 0 was dehned in Sec. since otherwise 


the minimum oPE^ 
holds true. 


In 


61 


is attained by a ball 37 Theorem 3.2] and the statement of the theorem 


Theorems 5.1.1 and 5.1.5], it is proved that the functional E^ admits 
a minimizer u = xFr £ Aloo(ni), Ep C for any R > 0, and after a possible redehnition 
on a set of Lebesgue measure zero, the set Ep is regular, in the sense that it is a union of 
hnitely many connected components whose boundaries are of class C'^h/2 Ei,... ,Ef,i 
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C be the connected components of Fji. By Lemma 4.1 we have N < Nq, |Lfc| > (5o 
and diam Fj. < Dq for all 1 < /c < and for some A^o > 1 and some constants Dq, (5o > 0 
depending only on m. Furthermore, we have 


dist{Fi, Fj) > 2R for i / j, (4-10) 

since otherwise it would be energetically preferable to increase the distance between the 
components. In particular, if i? > Dq the family of sets Fi,..., Fjy C generates a 
generalized minimizer {ui ,..., un) of E^o by letting Ui := XFi- Indeed, we have 

N N N 

e(m) > inf E^{u) = ^E^{ui) = ^E^{ui) >^e{\Fi\) > e{m), (4.11) 

1-^1=™ “T “T “T 


and so all the inequalities in (4.11) are in fact equalities. Since E^oixFi) > for each 

1 < i < A^, from (4.11) we obtain that each set Ei is a minimizer of E^o in Aoo{\Ei\). By 
Lemma 14.41 each set F) is bounded and connected, and dEi are of class C°°. □ 


The arguments in the proof of the previous theorem in fact show the following relation 
between minimizers of the truncated energy E^ and generalized minimizers of Eoo- 


Corollary 4.6 (Generalized minimizers as minimizers of the truncated problem). Let m > 
0 and R > d, let u ^ Aoo{m) be a minimizer of E^, and let u = where Ui are 

the charaeteristic functions of the connected components of the support of u. Then there 
exists a universal constant Ri > 0 such that if R > Ri, then (ui,... ,ui\f) is a generalized 
minimizer of E^a in Aooim). 


Proof. We choose Ri = max{i? 0 )-Dq}, where Ro and Do are as in Lemma 4.2 
have E^{xFo) = -^oo(xFo) for every connected component Fq of the minimizer. W 
same argument as the one used in the proof of Theorem 4.5 this yields the claim. 


We now provide some uniform estimates for generalized minimizers. 


Theorem 4.7 (Uniform estimates for generalized minimizers). There exist universal con¬ 
stants Jo > 0 and Dq > 0 such that, for any m > rho, where rho is defined in Sec. the 
support of each component of a generalized minimizer of E^o in Aooim) has volume bounded 
below by Jq o-nd diameter bounded above by Dq (after possibly modifying the components on 
sets of Lebesgue measure zero). Moreover, there are universal constants C, c > 0 such that 
the number N of the components satisfies 


cm < N < Cm. 


(4.12) 
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Proof. Let m > rho and let (xFd • • • j XFjv) t>e a generalized minimizer of E^o in Aaoim), 
taking all sets Fi to be regular. By |37[ Theorem 3.3] we know that there exists a universal 
iii -2 > such that 


< ^2 for alH G {1,..., N} . (4-13) 

Then by [37] Lemma 4.3] and the argument of |37[ Lemma 4.1] we have 

diam(Tj) < Dq, (4-14) 

for some universal Do > 0. On the other hand, we claim that taking R > Dq we have that 

N 

u{x) :='^XFi{x + AiRei), (4.15) 

where ei is the unit vector in the first coordinate direction, is a minimizer of in A(^{m). 
Indeed, since the connected components of the support of u are separated by distance 2R, 
we have 


N 


N 


= ^E^{xf,) = ^E^{xFi) = e{m) 


(4.16) 


2=1 


2=1 


At the same time, by the argument in the proof of Theorem 


4.5 


we have inf 


uGAocim) 


E^{u) = 


e{m) for all R sufficiently large depending on m. Hence, u is a minimizer of E^ in Aao{m) 


for large enough R. The universal lower bound \Fi\ > 6o then follows from Lemma 4.1 and 
our assumption on m. 


Finally, the lower bound in (4.12) is a consequence of (4.13), while the upper bound 
follows directly from the lower bound on the volume of the components just obtained. □ 


4.3 Properties of the function e(m) 

In this section, we discuss the properties of the functions e(m) = inf^^^ ^^^Eoo{u) and 
f{m) = e(m)/m, in particular their dependence on m. 

We start by showing that e{m) is locally Lipschitz continuous on (0, oo). 

Lemma 4.8 (Lipschitz continuity of e). The function e{m) is Lipschitz continuous on 
compact subsets of (0, oo). 

Proof. Let m, m! G [mo, mi] C (0, oo) and let (ui,... ,un) he a generalized minimizer of E^o 
in ^oo (m). For A = (m'/m)^/^, we define the rescaled functions with uf{x) = Ui{\ ^x). 


21 






For sufficiently large i? > 0, we define G Aoo{iti') by u^{x) := 
where ei is the unit vector in the first coordinate direction. We then have 


_ ^ r ^ r r 

E^{u^) = x^y \vui\dx + x^y / 

JR? JR 


Ui{x)ui{y) 
8vr|x - y\ 


dx dy + g{R), 


(4.17) 


where the term g{R) refers to the interaction energy between different components n^, , 

i A h of Clearly, we have g{R) —)• 0 for R —)• oo. It follows that 

N ^ N 


Ui{x)ui{y) 


Eooy)-e{m) <y - l\'y [ \Vui\dx+\X^ -l\y [ / „ I I 

y^JR^ y^jR^jR^ 87r|x-?/| 


dx dy+g{R). 

(4.18) 


In the limit R —?■ oo, this yields e{m') < Eoo{u^) < e(m)(l + Clm — m'\) for a constant 
C > 0 that depends only on mo, mi. Since m,m' are arbitrary and since e(m) is bounded 
above by the energy of a ball of mass mi, it follows that e is Lipschitz continuous on 
[mo,mi] for all 0 < uiq < mi. □ 

We next establish a compactness result for generalized minimizers. 

Lemma 4.9 (Compactness for generalized minimizers). Let m^ be a sequence of positive 
numbers converging to some m > mo, where mo is defined in Sec. as fc —>■ oo, and let 
{uk,i, ■ ■ ■ ,Uk^Nk) be a sequence of generalized minimizers of Eoo in Aoo{mk). Then, up to 
extracting a subsequence we have that Nk = iV G N for all k, and after suitable translations 
Uk,i Ui in BV (M^) as k ^ oo for all i G {1,..., N}, where (ui,... ,Un) G Aoo{m) is a 
generalized minimizer of E^ in Aoo{nn). 


Proof. By Theorem 4.7, we know that Nk < M G N for all k large enough. Hence, upon 
extraction of a subsequence we can asume that Nk = N for all k, for some N gN. For any 
f G {1,..., A^}, we also have 


sup 
k . 


\Vukfi dx < sup EooixB 

m£X " 


,1/S 


) < oo. 


(4.19) 


Moreover, again by Theorem 4.7 we have mk,i > ^o and supp(uA;^j) C H£)g(0), after suitable 
translations. Hence, up to extracting a further subsequence, there exist m^ > 6q and 
Ui G Aao(mi) such that mk,i —t m* and Uk,i Ui in BV(R.^), as fc —)• oo. Passing to the 
limit in the equalities m^ = X^i^i uik,i and e{mk) = obtain that 


N 

m = y^ mi 
i=l 


N 


and 


e(m) = ye{mi), 


(4.20) 


2=1 


where we used Lemma 4.8 to establish the last equality. Finally, again by Lemma [4.8| and 
by lower semicontinuity of E^q we have e{mi) < E^{ui) < liminffc_,.oo e(mfc,i) = e{mi), 
which yields the conclusion. □ 
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With the two lemmas above, we are now in a position to prove the main resnlt of this 
snbsection. 


Lemma 4.10. The set I defined in (2.17) is eompact. 


Proof. Since I is bonnded by |37[ Theorem 3.3], it is enongh to prove that it is closed. Let 
TTifc —7- m > 0, with nik G T, and let G ^oo(^fc) be snch that Eoo{uk) = e(mfc) for all 
A; G N, i.e., let Uk be a minimi zer o f the whole space problem with mass mk- We need to 
prove that m € I. By Lemma 4.9 there exists a minimizer u G Aoa{m) snch that Uk ^ u 
weakly in and Uk ^ u strongly in L^(]R^). In particnlar, there holds Eoo{u) = e(m) 

and hence m G X. □ 


Finally, we establish a few farther properties of e{m). 

Lemma 4.11. Let and A“ be the supremum and the infimum, respectively, of the 
Lagrange multipliers in ( |4.9[ ), among all generalized minimizers of Eoo with mass m > 0. 
Then the function e{m) has left and right derivatives at each m G (0, oo), and 


lim ~ = A- < A+ = lim ~ f ~ , 

h^0+ h m - m ^ 


(4.21) 


In particular, e is a.e. differentiable and e'{m) = A^ = A+ =: Am for a.e. m > 0. 


Proof. First of all, note that for m < mo, where mg is defined in Sec. the fnnction 
e(m) = m/(m) is given via (2.19), and the statement of the lemma can be verified explicitly. 
On the other hand, by definition we have A“ < A^. Fix m > friQ and let (ui,... ,un), with 
Ui = XFi, be a generalized minimizer of i?oo with mass m. We first show that 


Am > lim sup 
fe—>■0+ 


e(m + h) — e(m) 
h 


and 


. 4 - , ^ e(m) — e(m — h) 

Xt < lim inf 


h,—>0+ 


Indeed, for h > 0 let = Xf^ with so that \F^^\ = {^)\Fi\ 

(^)^/3 = l + 3 ^ + o(h), wehave 


(4.22) 

Since 


E^{uf) = EM + ^ 

h 


’dFi 


k{x) (x • f{x)) dP?{x) 


+ 


'■ f {x ■ f{x)) 

Unm Jqf^ Jf, \x - y\ 


dy dTL^iyx) + o{h), 


(4.23) 


where f{x) is the outward unit normal to dFi at point x. In view of the Euler-Lagrange 
equation (|4.9|), we hence obtain 


(4.24) 


\ 7 p 

Eoo{v!1) - EM = ^ (x • Hx)) dn‘^{x) = A ( \Ffi\ - \Fi\ ) + o{h), 
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where A is the Langrage multiplier in (4.9). Passing to the limit as h ^ 0"*“, this gives 

N 


lim sup 

/x—>-0+ 


e(m + h) — e{m) 
h 


< limsup-(^^.Eoo(M^) - < 

i=l i=i 


A. 


(4.25) 


Since (4.25) holds for all generalized minimizers, this yields the first inequality in (4.22). 
Following the same argument with h replaced by —h, and taking the limit as h —)• O"*", we 
obtain the second inequality in (4.22). 


Now, by Lemma 4.8 the function e{m) is a.e. differentiable on (0, oo), and at the points 
of differentiability we have e'(m) = A“ = A^ =: Am- Hence, for any h > 0 there exists 
m/j G (m, m + h) such that e is differentiable at rrih and 


e{m + h) — e{m) 
h 


> e'{mh) = Xm^ , 


so that 


lim inf 

/i^0+ 


e{m + h) — e{m) 

T 


> A := lim inf A. 

/i—>^0+ 


rrih ■ 


(4.26) 


(4.27) 


Let hk —>■ 0^ be a sequence such that Xm^^ —)• A as A; —>■ oo. If {u \,..., n^) are generalized 
minimizers with mass mh^. then by Lemma 4.9 they converge, up to a subsequence, to a 
generalized minimizer with mass m. In view of Lemma 4.4, up to another subsequence 
we also have that the boundaries of the components of the generalized minimizers with 
mass m/ij, converge strongly in to those of the limit generalized minimizer with mass 
m. Therefore, by (4.9) we have that A is the Lagrange multiplier associated with the limit 
minimizer. It then follows that A > A“, so that recalling (4.22) and (4.27) we get 


eim + h) — e(m) 
hm -^-= A^ 


h-5-0+ 


h 


(4.28) 


This is the first equality in (4.21). The last equality in (4.21) follows analogously by taking 
the limit from the other side. □ 


Remark 4.12. From the proof of Lemma 4.11 it follows that A^ are in fact the maximum 


and the minimum (not only the supremum and the infimum) of the Lagrange multipliers 
in (4.9), i.e., that A^ are attained by some generalized minimizers with mass m. 

Corollary 4.13. The function e{m) is Lipschitz continuous on [mo,oo) for any niQ > 0. 

Proof. This follows from ( 4.21| ), noticing that for all m > mo there holds 


—oo < 


inf .A-, <A-<A+< 




S sup 




(4.29) 


where M > 0 is such that I C [0, M], and we used (4.9) together with the uniform regularity 
from Lemma 4.4 for the components of the generalized minimizers. □ 
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4.4 Proof of Theorem 13.21 

In lieu of a complete characterization of the function f{m) and the set X, we show that 
f{m) is continuous and attains its infimum on X. 


The next result follows directly from Theorem 4.5 Theorem 4.7 and [37| Theorem 3.2]. 


Lemma 4.14. There exists a universal constant 5o > 0 such that for any m G (0, oo) there 
exist N > 1 and mi,... , G X such that mi > min{(5o, m} for all i = 1,..., N and 


N 

m = mi 
i=l 


N 


and 


f{m) = ^ —f{mi 


(4.30) 


i=l 


Theorem 3.2 is a corollary of the following result. 


Theorem 4.15. The function f{m) is Lipschitz continuous on [mo,oo) for any mo > 0. 
Furthermore, f{m) attains its minimum, i.e.. 


X* := < m* G X : f{m*) = inf f{m) 7 ^ 0. 

mSl 


Furthermore, we have f{m) > f* for all m > 0 and 


lim f{m) = 00 , 


m—^-O 


lim /(m) = f* 

m^oo 


lim ll/'l 


(m,oo) 


= 0 . 


(4.31) 


(4.32) 


Proof. Since f{m) = e{m)jm, the Lipschitz continuity of f{m) follows from Corollary 


4.13 By the continuity of /(m) and since X is compact, it then follows that there exists a 
(possibly non-unique) minimizer m* > 0 of f{m) over X. 


Turning to (4.32), the first statement there follows from (2.19). Let now u* = XF* £ 


•AooiiFi*) be a minimizer of E^o with m = m* for some m* G X*. Given fc G N, we can 
consider k copies of F* at sufficiently large distance as a test configuration. We hence get 
f{km*) < f{m*) for any A: G N, which implies f{m*) > lirninfm-^-oo/(™.). On the other 


hand, since f{m*) < f{m) for all m G X, by Lemma 4.14 we obtain 

N 


2=1 


/(m) = V —/(mi) > V —/(m*) =/(m*) Vm>0, 

^ ^ 777, ^ ^ m 


N 


2=1 


m. 

m 


(4.33) 


which gives the second identity in ( |4.32 ). Finally, by Corollary 4.13 we have 

e'{m)m — e(m 


lim |/^(m)| = lim 


m^ 


,. f {'Fl ) + 2||e c«) ^ 

< lim -^^ = 0 , 


m 


which yields the third identity in (4.32) 


(4.34) 

□ 
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5 Proof of Theorems 13.31 and 13.5 

5.1 Compactness and lower bound 


In this section, we present the proof of the lower bonnd part of the F-limit in Theorem 3.3 


Proposition 5.1 (Compactness and lower bonnd). Let {ue) G As, let be given by (2.4), 
let Vs be given by (|2.7|) and suppose that 


limsupe (us) < oo. 

£—^0 


(5.1) 


Then the following holds: 


i) There exists /i G A4^(T) n TT and v £ TL sueh that upon extraction of subsequences 
we have ^ /x m AI(T) and ^ v in TL. Furthermore, 


—Av = n — X 

a) The limit measure satisfies 


in V'{T). 


(5.2) 


Eoiia) < liminfe 

£-S>0 


(5.3) 


Proof. The proof proceeds via a seqnence of 4 steps. 


Step 1: Compactness. Since J^dp, 
and a snbseqnence snch that ps 
bonnd 


e = A, it follows that there is /x G At"*“(¥) with fj.dp = X 
p in A4(T). Furthermore, from (5.1) we have the uniform 


1 

2 


iVfgp dx 



y) dpe{x) dpeiy) < S ^/^EeiUe) < C. 


(5.4) 


By the definition of the potential, we also have fj. dx = 0. Upon extraction of a further 
subsequence, we hence get Vg ^ v in TL. Since ps ^ p in A4(T) and since the convolution 
of G with a continuous function is again continuous, we also have 


[ ( [ G{x-y)<p{x)dx') dpeiy) ^ [ ( [ G{x - y)(p{x) dx') dp{y) Vy? G P(T). (5.5) 
Jt \Jt J jt \jt J 

This yields by Fubini-Tonelli theorem and uniqueness of the distributional limit that 


v{x) 


G{x - y) dp{y) 


Furthermore, since satisfies 


-AVe = Pe - X 


for a.e. x G T. (5.6) 


in V'{T), (5.7) 
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taking the distributional limit, it follows that v satishes (5.2). In particular, (5.2) implies 
that /i dehnes a bounded functional on Ti, i.e. /r G . 


Step 2: Decomposition of the energy into near field and far field contributions. We split the 
nonlocal interaction into a far-held and a near-held component. For p G (0,1) and x G T, 
let rip{x) := ri(\x\/p), where p G C°°(R) is a monotonically increasing function such that 
p{t) = 0 for t < ^ and p{t) = 1 for t > 1. The far-held part Gp and the near-held part Hp 
of the kernel G are then given by 


Gp{x) = Pp{x)G{x), 


Hp-.= G- Gp. 


(5.8) 


For any u G ^e, we decompose the energy accordingly as Ee = + E^\ where 

g-4/3j^T)(^) _ 1-4/3 f f Q— y)u{x)u{y) dx dy 


T JT 


'^De^\u) = £ / |Vu| dx-|--e / / Hp{x — y)u{x)u{y) dx dy 


(5.9) 


T JT 


In the rescaled variables, the far held part E^^'^ of the energy can also be expressed as 


= 1- [ [ Gp{x-y)dpeix)dpe{y), 


(5.10) 


T JT 


where ps is given by (2.4). For the near held part E^'^ of the energy, we set := e 
and dehne h : —)■ M by 


ii(x) := u[xlte), (5-11) 

where is a torus with sidelength (cf. Sec. |^. In the rescaled variables, we get 

^-4/3^p)(^) _ ^1/3 j f f f e^^^Hp{£^^^{x — y))u{x)u{y) dxdy\ . 

2 J 

(5.12) 


Step 3: Passage to the limit: the near field part. Our strategy for the proof of the lower 
bound for (5.12) is to compare with the whole space energy treated in Section]^ and 
use the results of this section. We claim that 


liminfe > (1 — cp)Xf*, 

£—^0 

for some universal constant c > 0. 


(5.13) 
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Let r(x) := , X G be the Newtonian potential in and let r^(x) := , 

X G T, be the restriction of r(x) to the nnit torns. We also dehne the corresponding 
trnncated Newtonian potential : T —>■ M by 

T*{x) := (1 - r]p{x))T*{x). (5.14) 

By a standard resnlt, we have 

G{x) = T^{x) + R{x), X G T, (5.15) 


for some R G Lip(T). Hence 

Hp{x) = (1 - r]p{x))G{x) > (1 - r/p(x))(r#(x) - ||i?|Ucx>(T)) 

> (1 - 77 p(x))r#(x)(l - 47 r/ 3 ||i?||ioo(TP)) = (1 - cp)T*{x), (5.16) 

where c = 47r||i?||^oo(T)- Inserting this estimate into (5.12[), for cp < 1 we arrive at 


-4/3^p)(^ 

1 — cp 


>.1/3 




|Vtt| dx + ^ f I — y))u{x)u{y) dx dy 


= eV3 


(f \Vu\dx+ f [ — — u{x)u{y) dx dy\ . 

87 r|x-y| J 

(5.17) 


Next we want to pass to a whole space sitnation by extending the fnnction u periodically 
to the whole of and then trnncating it by zero ontside one period. We claim that after a 
snitable translation there is no concentration of the periodic extension of u, still denoted by 
u for simplicity, on the bonndary of a cnbe Qi^ := { — More precisely, we claim 

that 


'dQe, 


u{x — X*) d'H^{x) < 6 A, 


(5.18) 


for some x* G Qi^. Indeed, by Fnbini’s theorem we have 


A4 = 


/ udx= T-L‘^{{u{x) = 1} D {x ■ ei = t}) dt, 

JQt. J-ie. 


— If 


(5.19) 


where ei is the nnit vector in the hrst coordinate direction. This yields existence of Xi G 
(—snch that ?7^({n(x) = 1} H {x • ei = x)^}) < A. Repeating this argnment in 
the other two coordinate directions and taking advantage of periodicity of u, we obtain 
existence of x* G Q 4 snch that (5.18) holds. 
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Now we set 


u{x) := 


u{x — X*) 

0 


X G 
X G 


(5.20) 


We also introduce the truncated Newtonian potential on by 

1 - r]p{x) 


Tp{x) := 


47r|a:| 


X G 


(5.21) 


By (5.18), the additional interfacial energy due to the extension (5.20) is controlled: 
/ \Vu\dx= / \Vu\dx— / udx> / \Vu\dx — 6X. 


(5.22) 




'dQe, 


We hence get from (5.17): 


1 — cp 
> 


>ei/3 

A 

Imaudx 


' \Vu\dx + - 

R3 2 

|Vm| dx + - [ 


r^-i/ 3 „(x - y)u{x)u{y) dxdy - 6 A 


Tpg{x — y)u{x)u{y) dx dy\ — GAe^^^, (5.23) 


for any po > 0) provided that e is sufficiently small (depending on pq). By Corollary 4.6 


and Theorem 4.15 the first term on the right hand side is bounded below by Xf* as soon 
as po ^ .^ 1 - Therefore, passing to the limit as e —)• 0, we obtain (5.13). 

Step 4- Passage to the limit: the far field part. Passing to the limit ^ p in At(T), for 
the far field part of the energy we obtain 


linie \ [ [ Gp{x - y) dp{x) dp{y). 


£—^0 


(5.24) 


T JT 


At the same time, by (A.13) in Lemma A.2 in the appendix the set {{x,y) G T : a: = y} is 
negligible with respect to the product measure p0p on TxT. Therefore, since Gp{x — y) 
G{x — y) as p —)• 0 for all x 7 ^ y, by the monotone convergence theorem the right-hand 
side of ( |5.24 ) converges to fj fjG(x — y) dp{x) dp{y). Finally, the lower bound in (5.3) is 
recovered by combining this result with the limit of(5.13)asp—>-0. □ 


5.2 Upper bound construction 

We next give the proof of the upper bound in Theorem |3.3[ 

Proposition 5.2 (Upper bound construction). For any p G A4^(T) {^FL' with J^dp = X, 
there exists a sequence (n^) G Ae such that 

Pe ^ p in M.{T) and Vg ^ v inPL, (5.25) 
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as e —>■ 0, where fig, Vg and v are defined in (2.4), (2.7) and (3.4), respeetively, and 


limsupe ^^^Eg{ug) < Eo{fa). (5.26) 

€—)'0 

Proof. We first note that the limit energy is continuos with respect to convolutions. In 
particular, we may assume without loss of generality that dfj,{x) = g{x)dx for some g G 
C°°(T), and that there exist C > c > 0 such that 

c < gix) < C for all x G T. (5.27) 


We proceed now to the construction of the recovery sequence. For h > 0, we partition T 



Ug,s(x} ■■= u* for x G T, (5.28) 

i=i 

as the sum of Ng^^ suitably rescaled minimizers of Eoo{u)/ J^sudx. Note that 
u*(e~^^^x} dx = em*. To decide on the placement of we denote the number 

^ (i) 

of the centers in each cube as i.e., 

N^'i:=if{je{l,...,Ng,s} : G Qf}. (5.29) 


With this notation we have Ng^s = provided that supp(u£^, 5 ) H dQf = 0 for all 

i. The measure g is then locally approximated in every cube Qf by “droplets” uniformly 
distributed throughout each cube. Namely, we set 


and choose a 


U) 

e,d 


SO that 


pjQf) 1 


< dg, < 


(5.30) 


(5.31) 


(i) 

where dg^s ■= rniuj^j \afig — ay is the minimal distance between the centers, for some 
K' > K > 0 depending only on g. We also set 





if ag G Qt 


(5.32) 
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Then, if e is sufficiently small depending only on 5 and we find that G As for £ 
sufficiently small depending only on 5 and //. 

Finally, we define the measure ^^,5 associated with the test function constructed 
above, dfj,s, 5 ix) ■= £~‘^^^UE^six) dx, as in (2.4) and choose a sequence of <5 —>■ 0. Choosing a 
suitable sequence of e = —>■ 0, we have fJ-es,s A* in Af(T). For simplicity of notation, 
in the following we will suppress the (5-dependence, e.g., we will simply write instead of 
^£ 5 ,( 5 ) etc. 


It remains to prove (5.26). As in the proof of the lower bound, for a given p G (0,1) we 
split the kernel G into the far field part Gp and the near field part Hp. Decomposing the 
energy into the two parts in (5.9) and using (5.10), we have 


£ = ^ [ [ Gp{x-y)dpe{x)dpe{y)- 


(5.33) 


T Jf 


Since /ig ^ /x in A4(T), we can pass to the limit e —)• 0 in (5.33). Then, since the limit mea¬ 
sure p, belongs to Td', by the monotone convergence theorem we recover the full Coulombic 
part of the limit energy Eq in (3.1) in the limit p —)• 0. 

For the estimate of the near field part of the energy, we observe that 

Hp{x) < (1 -h cp)rf{x), (5.34) 

for some universal c > 0 (cf. the estimates in ( 5.16[ )). With this estimate, we get 

£-4/3^|2)(^^) < ^-1/3 |Vu^| dx-F ^e“'^/^(l-F cp) [ [ T*{x - y)ue{x)ue{y) dx dy 

Ue{x)Ue{y) 


T JT 


<£-^/^ \VUe\dx + £-^/\l + Cp) , , ,, , 

Jt JtJbi. (x) 87r|x-p| 


-he ^/^{l + cp) 


Ue{x)Ue{y) 


tJBp(x)\Bx {x) 87r|a:-y| 


dy dx 

dydx. (5.35) 


(i) 

By the optimality of u* and the fact that all 6^ ^ > 1, we hence get 


Jt JtJbip {x) 87r|a:-y| 


dydx < (1 -h cp)(A -h 0^(1))/*, 

(5.36) 


where the 0 ^( 1 ) term can be made to vanish in the limit by choosing £s small enough for 

(j) 

each 6 to ensure that all 9^ ^ —>■ 1. Since we can choose p > 0 arbitrary, this recovers the 
first term in the limit energy Eq in (3.1). 
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It hence remains to estimate the last term in (5.35). We hrst note that 


Ue{x)Ue{y) 


e / / "t v /-cxaj ^ 2 ; < A sup 


T JBp(x)\Bi (x) 


k - y\ 


dl^ejy) 

x£TJBp{x)\Bi (x) I® y\ 


(5.37) 


To control the last term, for any given x G T we introduce a family of dyadic balls Bj. ■ = 
B 2 -kp{x), k = 0,1,.... By ( |5.31[ ), we have Bp{x)\Bi^^{x) C \JklQBk\Bk+i for := 


\log 2 {p/de)~\ < 1 + log 2 (p/ci£), or, equivalently, 2~^^p > provided that e is sufficiently 
small depending only on 6 and p. Therefore, with our construction we have ^^(Bk) < 
2~^^Cp^ for some C > 0 depending only on p and all 0 < A; < K^. This yields 


sup / 

x&tJBp{x)\Bi. (x) 


dpejy) 

\x - y\ 


< 


Kp 

E 

fc =0 ' 


dpejy) 

,\Sfe+i \^-y\ 

- ^ p 


k=0 ’ k=0 

Since we can choose p > 0 arbitrarily small, this concludes the proof. 


2Cp^ ^ 8Cp^ 
4fc - 3 


(5.38) 


□ 


Remark 5.3. We note that the construction in Proposition |5.2| still yields, upon extraction 
of a subsequence, a recovery sequence for a given sequence of e = £„ 0. 


5.3 Equidistribution of energy 

We now prove Theorem |3.5| First, we observe that 


1 


dVe=£ ^/^\VUe\dx + -Vedpe, 


(5.39) 


where p^ is dehned in (2.4). We claim that the following lower bound for measures Ve-, given 


X G T and 5 G (0,1), holds true: 


liminl UeiBs{x)) > \Bsix)\Xf* 
£->0 


(5.40) 


As in (5.8), we split G into the far held part Gp and the near held part Hp, for some hxed 


p G (0,5). Since supp(Rp) C Bs{0), we obtain 

v^{Bs{x)) = ( \Vue\dx + ]-e~'^/^ [ [ Hp{x - y)ue{x)ue{y) dy dx 

JBs{x) ^ JBs{x) JBs{x) 

+ h / / Gp{x - y)dpe{y)dpe{x). (5.41) 

^ Jbs{x) Jt 
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Then, since Gp is smooth and HeiT) = A, by Corollary 3.4 the integral Jijp Gp{x — y) dfi^{y) 


converges to A fj. Gp{y) dy nniformly in x G T as e —?• 0. At the same time, by the definition 
of G and ( |5.34D we have 0 = Jj G{y) dy = fj Gp{y) dy + Hp{y) dy < fj Gp{y) dy + Gp^ 
for some nniversal C > 0. Hence, we get 


i>^{Bs{x))>e [ \Vue\dx + l- 

JBs(x) ^ 


2^ 


-4/3 


>Bs{x) JBs{x) 


Hp{x - y)ue{x)ue{y) dy dx - CXp^, 

(5.42) 


for £ snfficiently small and C > 0 nniversal. 

We now identify with its periodic extension to the whole of By Fnbini’s theorem, 
for a given 5' G (0, (5), there is t G {5', J) snch that 


[ ueix)dn^{x)<j^ [ U 

IdBtix) ^ - 0 Js' \JdBs{x) 


Ue(x) dT-L^iyx) ds = 


1 


5-6' 


Up dx. 


' Bs(x)\Bs,{x) 


(5.43) 


We then define Up G i?H(M^; {0,1}) by Up = UpXBt(x)- Recalling again Corollary 
obtain 


3.4 


we 


\\lup\dx= / \Vup\dx+ / Up{x)dld?{x) 
J Bt{x) JdBt{x) 


< 


[ \Vup\dx + GX5^e^/^, 
Jbs(x) 


(5.44) 


for some nniversal C > 0, provided that e is snfficiently small. We note that Up{x) < Up{x) 
for every x G Fnrthermore, for snfficiently small 6 we have Hp > 0 and 


Hp{x - y) > (1 - cp)V{x - y) for all |x - y| < \p, 


(5.45) 


for some nniversal c > 0 (where F is the Newtonian potential in as above). From (5.42|), 
(5.44) and (5.45) we then get 


i^eiBsix)) > e-^/3 / |Vn,| dx + L^e-4/3 

./to3 2 


'r3 JBp/ 2 {x) 


r(x - y)up{x)up{y) dy dx - CXp^, 

(5.46) 


for e small enongh. Letting now Up{x) := Up^e^^^x) be the rescaled fnnction which satisfies 
f Updx = - f Up dx = X\Bt{x)\£~^^^ + o(£~^^^), (5-47) 

Jr3 £ JrS 
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for every fixed po > 0 and e sufficiently small, we get 


Ve{Bs{x)) 

> (1 - cp)e^/^ 

(1 - 2cp)X\Bt{x)\ 


> 


/r3 Ue dx 


\Vue\dx + \[ [ T{x-y)ueix)ue{y)dydx] -CXp'^ 

Jm. 3 JBp^ix) J 

\Vue\dx + l- I I Tpgix-y)ueix)ue{y)dydx] -CXp‘^, 

JR3 J]R3 / 

(5.48) 


where Tpu is defined via (5.21). Recalling Corollary |4.6| and choosing po ^ -Rij we obtain 


liminf i^eiBsix)) > (1 - 2cp)Xf*\Bt{x)\ - CXp^, 

e-s>0 


(5.49) 


which gives (5.40) by first letting p —)• 0 and then 8' —)• 8. 


We now prove a matching upper bound. Notice that by the definition we have Ve{x) > 
C := — A| miupgT G(y)| for every x G T. Therefore, the negative part v~ of obeys 
^o^Xedps < ^\C\pe{U) for every open set U C T. In turn, since 
Vei^) = + 0£(1) by (3.8), it follows that the positive part of obeys i'^{U) = 

/f/nK>o} - ^f* + l\C\X + Oe{l). Hence l^sl = vf + ^7 is 

uniformly bounded as e —)• 0, and up to a subsequence I'e ^ v for some i' G A4(T) with 


z^(T) = Xf*. Since from the lower bound (5.40) we have I'iU) > A/*|17|, it then follows 


that du = Xf*dx. Finally, in view of the uniqueness of the limit measure, the result holds 
for the original sequence of e —)■ 0. □ 


6 Uniform estimates for minimizers of the rescaled energy 


In this section, we establish uniform estimates for the minimizers of the rescaled problem 
associated with over Ai from (2.23) and (2.24), respectively. The main result is a uniform 
bound on the modulus of the potential, independently of the domain size 


Throughout this section, T C with \F\ = XI is always taken to be such that U(, = XF 
is a regular representative of a minimizer of Ei over for a given A > 0 (for simplicity 
of notation, we suppress the explicit dependence of T on £ throughout this section). The 
estimates below are obtained for families of minimizers {ui^ as !'„—)• oo and hold for all 
U > ^ 0 ) where io > 0 may depend on A and the choice of the family. For simplicity of 
notation, we indicate this by saying that an estimate holds for £ 3> 1. 


Following M|64 we recall the notion of (A, ro)-minimizer of the perimeter (for a different 
approach that leads to the same regularity results, see |65]). 


34 










Definition 6.1. Given A, tq > 0 we say that a set F G is a volume-constrained (A, tq)- 
minimizer if 

P{F) < P{F') + A|FAF'| VF' C T^, s.t. {FAF') C Br^ and \F'\ = |F| , (6.1) 

where P{F) denotes the perimeter of the set F, and denotes a generic ball of radius r 
contained in Tg. 

The following result is a consequence of the regularity theory for minimal surfaces with 
volume constraint (see for instance Chapters III-IV], |61[ Section 4]). 

Proposition 6.2. Let F C he a volume-constrained (A,ro)-minimizer, with 

|F| G — Tq). Then dF is of class and there exist universal constants 6 > 0 

and c > 0 such that for all xq € F we have 

\For\ Br{xo)\ > cr^ /or a//r < min fro, , (6.2) 


where Fq is the connected component of F such that xq G Fq. 

Let r(x) := , X G be the Newtonian potential in and let T^{x) := 

Ti, be the restriction of r(x) to T^. Letting 


47r|ai| 




X G T^, 


by (5.15) we have for all i > 1: 

G(,{x) = vf (x) + Ri{x) for all x G T^, 
with Ri G Lip(T£) satisfying 

C C 

\Re{x)\ < ~j and |Vi?£(x)| < ^ for all x G , 

with a universal C > 0. 

Let now 

vf{x) := / Ge{x - y) dy, x G T^, 

Jf 

be the potential associated with F. Notice that vp satisfies 


, X G 

(6.3) 

(6.4) 

(6.5) 

( 6 . 6 ) 


—Avp = xf — -jj 2 and / vpdx = 0. 


(6.7) 


In particular, by standard elliptic regularity vp G C^’“(T£) for any a G (0,1) j^, and vp 
is subharmonic outside F, so that the maximum of vp is attained in F. Moreover, we have 
the following a priori bounds for vp throughout the rest of this section, vp always refers to 
the potential associated with the minimizer F). 
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Lemma 6.3. There exists a universal eonstant C > 0 such that 


-C <VF< C(A£)2/3, 


( 6 . 8 ) 


for all i ^ 1. 


is defined in (2.7), 


Proof. First of all, observe that vp{x) = for e = where 

in which /ig is given by (2.4) with Ue{x) = Furthermore, by a rescaling we 

have that is a minimizer of over As- Therefore, to establish a lower bound for vp, it 
is sufficient to do so for Ve- 

Let Gp and Hp be as in ( |5.8[ ) (with the choice of rj hxed once and for all), and note that 
there exists a universal pQ> t) such that Hp > 0 for all p G (0, po) and, hence. 


Ve{x)> / Gp{x-y)dp.e{y)- 

Jf 


(6.9) 


At the same time, by Corollary 3.4 and the boundedness of |VGp| we have 


Gp{x — y) dps{y) Gp{y) dy uniformly in x G T, 


( 6 . 10 ) 


as e —7- 0. Notice that from the dehnition of G we have 0 = Jrip G{x) dx = fj Gp{x) dx + 
Hp(x) dx. Therefore, by (5.15) we get 


-Gp'^ < 


[ Gp{x) 
Jt 


dx < 0, 


( 6 . 11 ) 


for some universal C > 0 and all p G (0, po)- Choosing p = min{/?o, A we then obtain 

> —2G for all e > 0 sufficiently small. 


On the other hand, by (6.5) there exists a universal constant C > 0 such that 

dy 


vpix) < G I < G 

JF\x-y\ \ 

< G{2pR^ + R-^\F\) 


r dy ^ \F\Bpix)\ 

ibr{x) \x-y\ R 


( 6 . 12 ) 


for any £ > 1 and i? > 0. The claim then follows by choosing R = □ 

Remark 6.4. Let Aq > 0 and let A G (0,Ao). Since vp > Amin^jg-ir G(a^)) h is also possible 
to obtain a lower bound on vp which depends only on Aq, and not on the family of the 
minimizers, provided that i > Iq for some .^o > 0 depending only on A. In this case all the 
estimates of this section still hold, but with constants that depend on Aq. 

We next obtain a pointwise estimate of the gradient of fi? in terms of vp itself. 
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Lemma 6.5. There exists a universal eonstant C > 0 such that for every 1 we have 


|V?;f(x)| < - (?;ir(x) + C) , 


(6.13) 


for any x G T^. 

Proof. Without loss of generality we may assume that x = 0. Arguing as in the proof of 
Lemma 6.3 and with the same notation, we can write 

|Vnir(0)|< [ \VGiiy)\dy = £ [\VG{y)\xF{yi) dy = [\VG{y)\dy^iy) 

Jf jt jt 

< ei/3 f \VGp{y)\dyM + f \VHp{y)\dy^iy), (6.14) 

JT JT 


where we recalled that s = £ Using (5.15), we have 

1-1 U I 


\VHp{y)\<{l + cp)\y\ Hp{y) + G\y\ p XB,\B^,.,{y), 


(6.15) 


for some universal c, C > 0 and all p G (0,po)- Substituting this into (|6.14) and recalling 
( [2^ and (|5^, we obtain 

|VnF(0)| ^ f Hp{y)dpeiy)+£~^^^ [ \y\~^Hp{y) u^iy) dy 

1+cp Mfi,i/3(0) 45,1/3(0) 

+ Ge^/^p~^ [ \y\~'^dpeiy)+s^^^ [ \'^Gp{y)\dpeiy) (6.16) 

JBpio)\B^^2(o) Jt 


< 


[ Hpiy)dpeiy) + C'{l + £^^^p ^A) + e^/^ [ \VGp{y)\dpe{y), 

JT JT 


for some universal G, G' > 0. Since by Corollary 3.4 and the smoothness of Gp we have 
JtIvgj X — y) \ dpe{y) —^ A /it |VGp(y)| dy uniformly in x G T as e —>■ 0, it is possible to 
choose eo > 0 sufficiently small independently of x such that the last two terms in the 
right-hand side of (6.16) are bounded by a universal constant for all e < eo. Thus, for every 
p < l/(2c) and e < eO) with eo depending on p, we have 


-\Vvf{0)\ <vf{0) + G - / Gp{y)dpe{y), 
n JT 


(6.17) 


where we also took into account that ^/^(O) = Jj Gp{y) dp£{y) + fj Hp{y) dpe^y). Finally, 
using (6.10) and (6.11), we obtain 


|Vx/r(0)| < -VFi0)+C{l + Xp^), 


(6.18) 


for some universal C > 0 and all e < eo, possibly decreasing the value of eo. The proof is 
concluded by choosing p < □ 
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Corollary 6.6. Let 1 and let x ^ F be a global maximum of vp. Then 


3 1 

VF{y) > T vf{x) - -C for all y € Bi/q{x) , 


(6.19) 


where C is as in (6.13). Furthermore, if vf{x) dx < C'\Br\ for some xq G T^, ^ — E 

and C' > 0; then 


VF{y) < C* + ‘iC for all y G Br{xo ), 


( 6 . 20 ) 


Proof. Since vf G C^iTg), for any y G Biiq{x) there exists 9 G (0,1) such that with the 
help of (6.13) we have 

vf{x) - VF{y) = Vvf{0x + (1 - 9)y) ■ {x - y) 

< \\Vvf{9x + {1 - e)y)\ 

6 

< ^VF{6x + {l-e)y) + ^C 




( 6 . 21 ) 


Similarly, letting y be a global maximum oi vf in Br{xQ) and letting xi G B^^xq) be 
such that ^^^(xi) = \Br\~^ vf{x) dx, we may write 

VFiy) < VF{y) - vf{xi) + C 

< |Vxir(6»xi + (1 - 6)y)\ \y - xi| + C' 

< ^ vf{0xi + (1 - 0)y) + + C" 


< 2 “*■ 2^ “*■ ’ 


which completes the proof. 


( 6 . 22 ) 

□ 


The next lemma provides a basic estimate for the variation of the Coulombic energy 
under uniformly bounded perturbations. 

Lemma 6.7. There exists a universal constant C > 0 such that for any i > 1 and for any 
F' C Ti, with FAF' C Br{xo) for some xq G and r > 0, there holds 


Vf dx — / vf' dx 


IF' 


< {‘^\\vf\\l'=°{Ti) + Cr'^) 


(6.23) 
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Proof. By direct computation, we have 


vpdx — / vpi dx 


IF' 



Te JTi 


x)Gi(x - y)xF{y) - XF'{x)Gi{x - y)xF'{y)) dxdy 


iXFix) + XF'ix))Gi(x - y)ixFiy) - XF'{y))dxdy 


JTi 


< 2 


Xf{x)G(,{x - y){xF{y) - XF'{y)) dx dy 


ft Jfi 


+ 


{xf{x) - Xf'{x))G(,{x - y){xF{y) - XF'{y))dxdy 


< 2 


ne 


JTe 

VF{y){xF{y) -XF'{y))dy 


+ 2 



Tf J Br{y) 


Gi{x - y){xF{y) - XF'{y)) dx dy 


^ (2||^^f||l°o(t^) + Gr"^) \FAF'\ 


for some universal C > 0, where we used (6.4) and (6.5) in the last line. 


(6.24) 

□ 


Lemma 


6.7 


implies that minimizers of Ei are volume constrained (A, ro)-minimizers of 
the perimeter for rg = 1 and A = ||'Cf||l°°(t^) + G, with G > 0 universal. In particular, by 
we get A < C'(A£)^/^, provided that £ 3> 1. Therefore, from Proposition 


6.3 


Lemma 


obtain the following result. 


6.2 


we 


Proposition 6.8. There exist universal constants c > 0 and <5 > 0 such that for all i ^ 1 
and all xq € F there holds 


|Fo n 7?r(a:o)| > cr^ for all r < 


(A£)2/3 ’ 

where Fq is the connected component of F such that xg € Fq. 


(6.25) 


We now show that the potential vp is bounded in by a universal constant as 

f —)> oo. 

Theorem 6.9 (L°°-estimate on the potential). There exists a universal constant G > 0 
and a constant £g > 0 such that for all £ > io we have 


\xf\\l°°{Ti,) < G. 


(6.26) 


Proof. Observe first that by (6.8) we have vp > —C, for some universal constant G > 0 
and £^ 1. Therefore, letting V := maxvp(x), the thesis is equivalent to showing that 

xGTe 


V <G, 


(6.27) 
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for some universal C > 0 and large enough i. 

We first prove ( |6.27[ ) with the constant depending only on A. Partition into N 
cubes of sidelength L = with chosen to be the smallest integer such that 

L < min where c and 5 are as in (6.25). Note that with our choice of L 

we have N > 216A^£^/(c(5^). If £ is sufficiently large (depending on A), we also have that 
6{\£)~‘^/^ < \L < In particular, any ball of radius 5{\£) can be inscribed 

into a union of 27 adjacent cubes of the partition and stay at least distance 5{\£)~‘^^^ from 
the boundary of that union. Hence, by (6.25) and a counting argument we get that at 
least cubes do not intersect F, so that we can find disjoint balls Hi,..., Bm of radius 
< g not intersecting F, with M > '^N. 

Recalling that vpdx = 0 and that vp is bounded below by —C, for ^ 1 we get 


0 = 


/ vpdx>^ [ 

^—1 Bi 


Vp dx — C£'^. 


(6.28) 


It follows that there exists an index i such that, for some universal (7^ > 0, we have 


vpdx < CM-^£^ < C'lHJ . 


(6.29) 


'Bi 


We then apply the second part of Corollary 6.6 with xq = Xi, where Xi is the center of Bi, 
to obtain 


|'yF(a;)| < C for all x G Bi, (6.30) 

for some universal (7 > 0. 

Let now x £ F he a global maximum of vp, so that vp{x) = V, and assume that 


n‘^{F n dBr{x)) > ^V\F n Br{x)\ for any r G (0, L/2). 


(6.31) 


Letting m(r) := \F n Br{x)\, so that = H^(F n dBr{x)) for a.e. r, (6.31) can be 


written as 


> ^Vm{r) for a.e. r G (0, L/2) 
dr 9 


(6.32) 


Integrating (6.32) over (ro,L/2), we get (for a similar argument, see the proof of 
Theorem 3.3]) 


m(ro) < m{L/2) 


Notice now that, as in Proposition |6.8| from Lemma 6.7 it follows that 

(5 


m{r) > cr'^ for all r < min I 1, 


V 


(6.33) 


(6.34) 
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In particular, if ro = bjV < L/4, we have 


^ < m(ro) < , 


(6.35) 


for some universal constant C > 0, which implies (6.27) with the constant depending only 
on A. 


On the other hand, if (6.31) does not hold, there exists r G such that 

n dBr{x)) < ^V\F n Br{x)\. 


(6.36) 


We claim that, as in the proof of Lemma 4.2, if (6.27) does not hold, it is convenient to 
move the set F 0 Br{x) inside the ball Bi. Indeed, we dehne L) := (x* — x) + (F H Br{x)) 
and u = U£ — XFnBr{x) + XFi- Note that by construction F f] Bi = 0, so u is admissible. 
By minimality of Ui and using (6.5), (6.30) and (6.36), we get 


Ei{ue) < Ei{u) 

= Ei{u£) + 2'H'^{F n dBr{x)) + / vpdx — 


'Fi 


' FnBr{x) 


vp dx 



GAx — y)dxdy+ / / GAx — y)dxdy 

I Fi JFr\Br{x) JFFBrix) JFnBrix) 


< EeiuA + {-V + G ] \E n Br{x)\- 


' Ff\Br{x) 


vp dx , 


(6.37) 


for some universal C > 0, provided that i ^ 1. Notice now that Corollary 6.6 implies that 


vp{x) > -V — G for any x G Br{x), 


(6.38) 


for a universal C > 0. Hence 


0 < -yj |FnH^(x)|, (6.39) 

for some universal C > 0 and i ^ 1, which leads to a contradiction if V is too large. 


Lastly, to establish (6.27) with C universal, we note that using (6.27) with the constant 
depending on A one gets that the density estimate in (6.25) holds for all r < tq with some 
ro > 0 depending only on A, for i ^ 1. We can then repeat the covering argument at the 
beginning of the proof with L > 0 universal, provided that 1, and obtain the conclusion 
by repeating the above argument. □ 

From Theorem 6.9 and the arguments leading to Propositionwe obtain an improved 
density estimate for minimizers of Ei. 
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Corollary 6.10. There exist a universal constant c > 0 and a constant > 0 such that 
for all xq € F and all i > io we have 

|Fo n Br{xo)\ > cr^ for all r < I, (6.40) 

where Fq is the connected component of F such that xq G Fq. 

Finally, we establish a uniform diameter bound for the connected components of the 
minimizers in Theorem 16.91 

Lemma 6.11 (Diameter bound). Let Fq be a connected component of F. Then there exists 
a universal constant C > 0 such that 


diamFo < C, 


(6.41) 


for all i ^ 1. 


Proof. Assume that diamFo > 2. Arguing as in the proof of Lemma 6.5 and using its 
notations, for any x G and a universally small po > 0 we have 

vf{x)> f — 1+ f Gp{e^/^x -y)dne{y), (6.42) 

for ah p G (0, po). Observe that by ( 6.10[ ) and (6.11) the last term in the right-hand side of 
(6.421 can be bounded below by —2C\p^, for £ S> 1 and C > 0 universal. Taking p < 
and using (6.26), we then get 


dy 


JF^Bnix) 1^; - y| 

with a universal C > 0, for any i? > 1 and x G T^, provided that £ 3> 1 independently of x. 


that 


Recalling (6.40) and arguing as in Lemma 

dy 


C > 


I FonBji{xo) 


|a^o - y\ 


< c, (6.43) 

r of X. 

for all 1 there exists xq G Fq such 
> c min{log (diam Fq) , log R} , (6.44) 


4.2 


for some universal c,C > 0. The claim then follows by choosing a universal R that is 
sufficiently large. □ 


7 


Proof of Theorem 




For A > 0, let (us) G As be a family of the regular representatives of minimizers of E^, and 
let Ne and £ i?R(M^;{0,1}) be as in the statement of the theorem. Without loss of 
generality we may set x^^k = 0 in the statements below. We need to show that there exists 
eo > 0 such that for ah e < £q: 
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i) There exist universal constants C,c> 0 such that 




/ dx > C£. 


(7.1) 


ii) There exist universal constants C,c> 0 such that 

supp(rie,fc) C Bq^i/ 3{0), < CXe~^^^. (7.2) 


hi) There exists a collection of indices such that (#/e)/A4 —?• 1 as e —)■ 0 and, upon 
extraction of a subsequence, for every sequence e„ —0 and every kn G Ie„ there holds 
hn —t h in L^(M^), where Un{x) ■= and u is a minimizer of E^o over 

Aaoitn*) for some m* € I*. 


The estimate for the potential in (i) follows from Theorem 6.9, setting G Ai^ 


with = £~^A and noting that with = XF we have vp = Ve{-lie)- Similarly, the volume 
estimate in (i) follows from Corollary 6.10 The inclusion in (ii) follows from Lemma 6.11 by 
a rescaling. The estimate for in (ii) follows from (i) and the fact that fjU^dx = 

We turn to the proof of statement (hi). Given h > 0, let > 0 be the number of the 
components such that for hg^fc(x) := ttg^fc(e^/^x), we have 


Eoo{Ue,k) > if* +d) U£,kdx. 


(7.3) 


By (3.8), (3.12) and the arguments in the proof of Proposition 5.1 we have, as e —)■ 0, 

Ne 

Xf* = £~'^/^Ee{Ue) +Og(l) > Eoo{Ue,k) + 0^(1) 


k=l 


> £ 


1/3 


+ ( 5 )^/ Ue,kdx + f* ^ / Ue,kdx) + Oe{l) 


k=l 


> Xf* + c SN^^s + 05 ( 1 ), 


(7.4) 


where we suitably ordered all Ug ^ and included a possibility that the range of summation 
is empty in either of the two sums. Hence, = o(e“^/^), and by (ii) it follows that 

Aig ,5 = o(A^g) for all <5 > 0. This implies that for every h > 0 there is £5 > 0 and a 
collection of indices leg satisfying {Hleg)/Neg —)• 1 such that Eoa{u£g^k)/ Ueg,k dx —)■ f* 
uniformly in k € I^g as h —)• 0. By (ii), for every sequence of 0 and every choice 

of kn G leg^ the sequence Un ■= li-eg ,k„ is supported in for some i? > 0 universal 

and equibounded in i?H(M^). Hence, upon extraction of a subsequence we have hn —)> u 
in L^(M^) with m := J^^udx > 0. At the same time, by lower semicontinuity of E^o we 
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4.15 


the latter is, in fact, an equality, and 
—)• 1, is a minimizing sequence for E^o 


also have Eoo{u)/m < f*. Then, by Theorem 
so Unix) := uiXnX) with An := (m“^ Ue^k dx)^ 
over Aooim) (cf. ( 4.17| )). Thus, ft is a minimizer of E^o over Aooim). Again, by Theorem 
MB we then have m G I*. □ 
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A Appendix 


We recall that by the Riesz-Fischer theorem, the space of signed Radon measures AI(T) is 
embedded in the space of distributions via the identification 





yip G C“(T). (A.l) 


On the other hand, any measure p G AI'^(T) C %' (recall 
extended by continuity to an element of the dual space Ti', 


the definition in (2.12)) can be 
which we still denote by /x, such 


that 


ifdp = Hiy>,y)i 


V(/p G ?7nC°(T). (A.2) 


Lemma A.l. Let p G AI'^(T) r\%' and u gT-L. Then, up to taking the precise representa¬ 
tive, u belongs to L^{T,dpi) and 


i{u,T)n'= / udfi. 


(A.3) 


Proof. The result follows as in 66 Theorem 1]. For the reader’s convenience we include a 


simple alternative proof here. Since u G TL,hy |67[ Section 4.8: Theorem 1] we can identify 
u with its precise representative and find a sequence Uk G 0 C'®(T) such that Uk ^ u in 
TL, and 


Ukix) 


uix) 


for all X ^ N, (A.4) 
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where C T is a set of zero inner capacity, that is, for any compact set K <Z N there 
exists a seqnence ipn G T-LCi C'^(T) snch that tpn ^ 0 in Ti and ipn = 1 on K. Since ^ £ Ti' 
we have ii{K) = 0 for all compact K G N, so that 


IJ.{N) = sup ^{K) = 0. 
KCN 


(A.5) 


Since the functions are continuous for all /c G N, we have 


n{uk,lj)H'= / Ukdn, 

JT 


Therefore, by (A.2) we get 

H (I Uk' I Oik^k' 1 ^k' I /^C^)) 

JT 


for all k' £N, where 


It then follows 


0^k,k' • / \'^k ■ 

JT 


(A.6) 


(A.7) 


(A.8) 


\uk - Uk'\\L^(T,dii) < lll^fc - Uk'\ - 0;k,k'\\n\\k-\\n' + M(’I’)||'Ufc - Uk'\\Li{T) 

= \\V{uk-Uk')\\L^T)\\lJ‘\\w + ^{'^)\\uk-Uk'\\L^(T). (A.9) 


Since Uk is a Cauchy sequence in Ti, hence also in L^(T), from (A.9) it follows that Uk is 
a Cauchy sequence in L^(T,d/i) and, therefore, converges to some u G L^{T,dfi). In fact. 


Therefore, from (A.6) we get 


passing to a subsequence and using (A.4) and (A.5), we have u{x) = u{x) for /x-a.e. x G T. 


H{u,k-)n'= lim = 1™ / Ukd^i= / ud^i 

k^oo k^oo Jj Jy 


which concludes the proof. 


(A.IO) 

□ 


The following lemma characterizes the measures in terms of the Coulombic potential, 
see |38[ Lemma 3.2] for a related result. 


Lemma A.2. Let ^ G M^{T)r\T-L', and let G :T ^ {—oo, Too] be the unique distributional 


solution of (2.6) with G(0) = +oo. Then the function 


v{x) := / G{x — y)d^{y) xGT 

JT 


(ATI) 
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belongs to % and solves 


— / vAipdx= ifdfi 

Jt Jt 


V(^ G c°°{T)nn. 


Moreover v G L^(T, d/i) and 


G{x — y) d^{x) d^{y) = / vdy,= / \Vvfdx. 

IT JT Jt Jt 


(A.12) 


(A.13) 


Proof. By the definition of G and the fact that G G L^(T), the function v belongs to L^(T), 


solves (A.12) and has zero average on T. On the other hand, by (2.12) one can define a 


functional G PL' such that T^{(p) = Jjipdy for every cp G (^“(T) nTi. Therefore, by 
Riesz Representation Theorem there exists v £PL such that 


— / vApdx = {p,v)'p = — vApdx 

Jt Jt 


yp G G°°iT)nn. 


(A.14) 


Thus, since A is a one-to-one map from G°°(T)riPL to itself, we conclude that v = v almost 
everywhere with respect to the Lebesgue measure on T and, hence, v £ PL. 

Let now p G C°°(T) be a radial symmetric-decreasing mollifier supported on Ri/8(0), 
let Pn{x) '■= n^p{nx), so that pn —t do hr D'(T), and let fn G (^“(T) be defined as 


fnix)-.= [ pn{x-y)dp{y) x G T. (A.15) 

JT 

Then, if the measures pn G Af'''(T) 0 PL' are such that dpn = fndx, we have 

in PL' and ^ /r in A4(T). Letting also Vn{x) := fjG(x — y)dpniy), we observe that 

Vn ^ V £ PL, and pn® Pn ^ P® p Af (T x T). For all M > 0, we then get 


Gm{x- y)dp{x)dp{y) = lim / / Gm{x - y) dpn{x) dpn{y) , (A.16) 

It Jt Jt Jt 

where we set Gm{x) '■= min(G(x),M) G C'(T). By Monotone Convergence Theorem we 
also have 


G{x - y) dp{x) dp{y) = lira. / Gm{x - y) dp{x) dp{y) 
It Jt m^°°JtJt 


(A.17) 


Recalling (A. 16), it then follows 


G{x - y) dp{x) dp{y) = lim lim / / Gm{x - y) dpn{x) dpn{y) 

M—^oo n—>-oo _y.Jp Jy 


< lim 


72—^OO J t/T 


G{x-y) dpn{x)dpn{y) 


= lim / Vndpn= lim \\vn\\n = \\v\\n ■ (A.18) 

n^oo Jy n^oo 
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Together with the fact that G is bounded from below, by Fubini-Tonelli theorem this implies 
that V G L^{T,dfi), with ||u||j, i(ir,d^) < \\v\\^. 

It remains to prove (A.13). We reason as in |68[ Theorem 1.11] and pass to the limit, 
as n —>■ oo, in the equality 


/ VndHn= / \VVn\^dx, 

h Jt 


(A.19) 


which holds for all n G N. Notice that the right-hand side of (A.19) converges since Vn ^ v 
in so that 


lim / \Vvn\^dx= / 
n^oo Jj Jy 


(A.20) 


In order to pass to the limit in the left-hand side of (A.19), we write 


/ Vndnn= / / G{x - y) d^in{x)d^n{y) = / / Gnix - y) dy{x)dy{y), (A.21) 

It Jt Jt Jt Jt 


where we set 


Gn{x) 

Pn{x) 


G{x - y)pn{y) dy, 


Pn{x - y)pn{y) dy. 


JT 

We claim that there exists C > 0 such that 

\Gn{x)\<G{l + \G{x)\) 


(A.22) 
(A.23) 

(A.24) 


for all X G T. Indeed, we can write G 


T^ -|- i? as in (5.15). Letting 


^t{x)-= [ ^*ix - y)pniy)dy and Rn{x) := [ R{x - y)pn{y) dy , (A.25) 

JT JT 

we have that i?„ —?■ R uniformly as n —?• oo. Moreover, since T"^, T^ and pn are periodic 
when viewed as functions on rewriting the integrals as integrals over subsets of and 
applying Newton’s Theorem we get 


Tt{x) 

r#(x) 


47r|x| / T*{x - y)pn{y)dy =\x\ j 

JSi/4(o) Jbi/4(o) |3;-y| 

pn{y) 


dy 


+ \x 


L 


, , dy< Pn(v) dy = 1 

/4(0)\S|,|(0) \y\ 


Pn{y) dy 


'Spi(O) 

for all k| < 7 - 
4 

(A.26) 
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Since also 


r*{x) 
r#(x) 


\ = ATr\x\[ T*{x-y)pn{y)dy <C for all ^ < |x| < (A.27) 

J ^Bi/8(0) 4 2 


this proves (A.24). 


From the fact that Gn{x) —)• G{x) for all x G T, by (A.24) and the Dominated Conver¬ 
gence Theorem we get 


lim / VndyLn= hm / / Gn{x - y) dp{x)dyL{y) 

n^oo n^oo 


G{x — y) dp{x)dp{y) = / v dp. 

ITJT Jt 


From (A.19), (A.20) and (A.28) we obtain (A.13) 


(A.28) 


□ 


Lemma A.3. Let G he as in Lemma A.2 and let p G Af'''(T) satisfy (2.13). Then p G 

M+{T)nn'. 

Proof. Let cp G C'^(T) n H. Using the same notation and argnments as in the proof of 


Lemma A.2 with the help of Canchy-Schwarz ineqnality we obtain 


ipdp= lim / ipdpn= lim / Vip-Vvndx 


< 


n lim { / G{x - y) dpn{x) dpn{y) 

n^oo yjj Jj 


= \m\'H lim / / Gn{x - y)dp{x)dp{y) 
n^oo yjj Jj 


= WWH 


G{x - y) dp{x) dp{y) 


T JT 


which yields the ineqnality in (2.12). 


(A.29) 

□ 
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